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ABSTRACT 
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Context. Herschel provides crucial constraints on the IR SEDs of galaxies, allowing unprecedented accuracy on the dust mass esti- 
mates. However, these estimates rely on non-linear models and poorly-known optical properties. 
• Aims. In this paper, we perform detailed modelling of the Spitzer and Herschel observations of the LMC, in order to: ( i) systematically 

study the uncertainties and biases affecting dust mass estimates; and to (ii) explore the peculiar ISM properties of the LMC. 
Methiods. To achieve these goals, we have modelled the spatially resolved SEDs with two alternate grain compositions, to study 
Q the impact of different submillimetre opacities on the dust mass. We have rigorously propagated the observational errors (noise and 

^ calibration) through the entire fitting process, in order to derive consistent parameter uncertainties. 

Results. First, we show that using the integrated SED leads to underestimating the dust mass by ^ 50 % compared to the value ob- 
tained with sufficient spatial resolution, for the region we studied. This might be the case, in general, for unresolved galaxies. Second, 
' ' we show that Milky Way type grains produce higher gas-to-dust mass ratios than what seems possible according to the element abun- 

dances in the LMC. A spatial analysis shows that this dilemma is the result of an exceptional property: the grains of the LMC have on 
T-H average a larger intrinsic submm opacity (emissivity index /3 - 1.7 and opacity /^abs(160 fim) = 1.6 m- kg"') than those of the Galaxy. 

^ By studying the spatial distribution of the gas-to-dust mass ratio, we are able to constrain the fraction of unseen gas mass between 

^10, and ^ 100 % and show that it is not sufficient to explain the gas-to-dust mass ratio obtained with Milky Way type grains. Finally, 
\^ we confirm the detection of a 500 /jm extended emission excess with an average relative amplitude of ^ 15 %, varying up to 40 %. 

This excess anticorrelates well with the dust mass surface density. Although we do not know the origin of this excess, we show that it 
T-H is unlikely the result of very cold dust, or CMB fluctuations. 

Key words, dust - ISM: abundances ~ Galaxies: ISM, dwarf, starburst - Individual: LMC 

• • 1. Introduction probe interstellar matter in regions where no other counterpart is 

^ accessible. 
k>( The infrared (IR) spectral energy distribution (SED) is widely 

used to derive the global properties of a system, such as its Unfortunately, there are several fundamental systematic un- 
instantaneous star formation rate, its dust and eventually gas knowns inherent to dust modelling, which are questioning the 
masses, and the compactness of the star forming region. The ad- reliability of these diagnostics. First, the microscopic proper- 
vent of the Herschel Space Observatory has opened the most ties of the grains are still poorly known. In the Milky Way, 
important spectral window to perfect these diagnostics, by ob- the most complete and accurate models are constrained by ob- 
serving the far-IR to submillimeter (submm) wavelengths (60- servations of high latitude cirrus clouds: their IR emission, ul- 
600 pm). Indeed, this regime samples the peak and Rayleigh- traviolet (UV)-to-near-IR extinction, and elemental depletions 
Jeans wing of the dust emission. It consequently constrains the (Zubko et al., 2004; Draine & Li, 2007; Compiegne et al., 201 1). 
emission by grains in thermal equilibrium with the radiation The authors performing these models derive the size distiibution 
field, present in the different phases of the interstellar medium and abundance of the different grain species - silicates, carbon 
(ISM), including the coldest, most massive components (down grains (graphite or amorphous carbons) and polycyclic aromatic 
to dust temperatures of Tjust ^ 12 K). This spectral domain is hydrocarbons (PAH). Zubko et al. (2004) demonstrated an im- 
therefore crucial to derive accurate dust masses, and physical portant degeneracy by presenting complete fits of the same data 
conditions, and can be used as a powerful, unprecedented tool to set (Galactic emission, extinction and depletion), with five dif- 
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ferent dust compositions, alternating bare and coated grains, as 
well as crystalline and amorphous solids. Thus, the derived dust 
properties depend on the assumed chemical composition of each 
species. The UV to millimetre (mm) opacities are sensitive to 
the grain composition. They are derived from sparse constraints 
including astrophysical features, laboratory spectra of analogs of 
interstellar dust materials, and theoretical solid state physics {e.g. 
Weingartner & Draine, 2001; Draine, 2003b). Their universality 
is doubtful. 

The second major source of uncertainties concerns the 
macroscopic variations of these microscopic grain properties, as 
a function of the environment. These variations are numerous; 
some are speculative: ( i) PAHs are known to be destroyed in H ii 
regions (e.g. Madden et al., 2006); (ii) the variations in the Ry 
parameter of the extinction curve is interpreted as variations of 
the grain size distribution (Draine, 2003a; Fitzpatrick & Massa, 
2005); (Hi) coagulation occurs in dense regions [e.g. Stepnik 
et al., 2003; Berne et al., 2007); (iv) blast waves are responsible 
for grain fragmentation and erosion in the low-velocity phase 
(Jones et al., 1996) and destruction close to the remnant (Reach 
et al., 2002); (v) the dust abundances and properties are thought 
to evolve with the metallicity of the ISM (Galliano et al., 2003, 
2005, 2008a); (vi) a transition from amorphous to crystalline sil- 
icates is observed in protostellar objects (e.g. Hallenbeck et al., 
2000; Poteet et al., 2011). This list is not exhaustive. Moreover, 
when considering the SED of a given region, it is possible to 
confuse variations of the physical properties of the grains (e.g. 
their optical properties) with variations of their physical condi- 
tions (e.g. the starlight intensity to which they are exposed). This 
problem becomes even more intricate, when considering the in- 
tegrated SED of a galaxy. 

The various processes controlling the lifecycle of dust 
throughout the ISM are not known with enough precision to 
break these kinds of degeneracies. Even the origin of interstel- 
lar dust is uncertain. The contribution of supemovae (SN) and 
asymptotic giant branch (AGB) stars to the observed content of 
ISM dust, and the dust growth in interstellar cloud is still de- 
bated (e.g. Galliano et al., 2008a; Draine, 2009). Dust is believed 
to constantly evolve throughout the ISM, being photoprocessed, 
altered by cosmic rays, accreting atoms in dense regions, and be- 
ing shattered in shocks (Jones, 2004, for a review). We are com- 
pelled to find observational cases where there will be enough 
redundancy in the data to isolate one of these processes. This is 
the goal of this paper 

The present article scrutinizes the different methodologi- 
cal biases, as well as the fundamental physical processes af- 
fecting the dust mass estimate in galaxies. Our demonstration 
is performed on the Herschel and Spitzer observations of the 
Large Magellanic Cloud (LMC; d - 50 ± 2.5 kpc; Schaefer, 
2008). Due to its proximity, it is an ideal laboratory to study 
the variations of the far-IR properties, down to spatial scales of 
^ 10 pc (SPIREsoo^jm angular resolution of 36")- Moreover, it 
offers an environment containing massive star clusters, allow- 
ing us to study the impact of intense star formation on the sur- 
rounding ISM. Finally, its metallicity is moderately sub-solar, 
with (0/H)lmc - 0.5 X (O/H)0 and (C/H)lmc - 0.3 x (C/H)o 
(Pagel, 2003). The comparison of its dust properties with those 
of the Galaxy therefore provides insights on cosmic dust evolu- 
tion. 

In general, low-metallicity dwarf galaxies, a category to 
which the LMC belongs, have peculiar dust properties. They ex- 
hibit a deficit of PAH strength, that appears to be correlated with 



the metallicity of their ISM. The origin of this trend is still de- 
bated: ( i) PAHs could be more massively destroyed by permeat- 
ing hard radiation, in sub-solar ISM (e.g. Madden et al., 2006); 
( ii) the delayed injection of PAHs by AGB stars could explain 
their lower intrinsic abundance, in young systems (Galliano 
et al., 2008a); (Hi) the PAHs could form in molecular clouds, 
which have a lower filling factor, in low-metallicity environ- 
ments (e.g. Sandstrom et al., 2010). 

The IR SED of dwarf galaxies usually peaks at shorter 
wavelengths, indicating hotter equilibrium grains, on average. 
In addition, the mid-IR continuum is steeply rising, similarly to 
what is observed in Galactic compact H ii regions (e.g. Peeters 
et al., 2002). This peculiar mid-IR continuum, was modelled by 
Galliano et al. (2003, 2005) with an increase of the very small 
grain abundances, which could be the consequence of the high 
number density of shock waves. It is consistent with the pecu- 
liar shape of the extinction curve, in the Magellanic clouds: their 
lower 2175 A bump, and their steeper near-UV rise could be the 
result of an excess of small grains (Weingartner & Draine, 2001; 
Galliano et al., 2003, 2005). This typical continumm shape was 
reported as a "mid-IR excess" compared to the Galactic SED by 
Bernard et al. (2008), in the LMC. 

The dust-to-gas mass ratio increases with metallicity 
(Lisenfeld & Ferrara, 1998; Draine et al., 2007; Galliano et al., 
2008a; Engelbracht et al., 2008). At first order, it can be un- 
derstood since dust is made out of the metals synthesized by the 
various stellar populations. However, the detailed dependency of 
the dust-to-gas mass ratio with the metallicity remains unknown. 
In addition, the gas mass, especially the molecular phase, is un- 
certain, because the CO-to-H2 conversion factor is a strong un- 
known function of the environment - in particular, it is a func- 
tion of the metallicity. A given CO line intensity will translate in 
a larger molecular gas mass in a lower metallicity environments 
(e.g. Madden et al., 1997; Israel, 1997; Leroy et al., 2007, 2009, 
201 1). Moreover, we can expect the presence of large quantities 
of H2 in regions where no CO at all is detected (Madden, 2000). 
Bernard et al. (2008) unveiled a "far-IR excess", compared to the 
gas column density, in the LMC, which is likely the evidence of 
such a gas reservoir. 

Finally, at submm wavelengths, the SED of dwarf galax- 
ies differs significantly from solar metallicity systems. Galliano 
et al. (2003, 2005) reported in four blue compact dwarf galaxies 
an excess emission at 850 ^m (SCUBA) and 1 .2 mm (MAMBO). 
Such an excess was then confirmed in other similar systems 
(Dumke et al., 2004; Galametz et al., 2009, 2010; Bot et al., 
2010; Grossi et al., 2010). It extends up to centimetre (cm) 
wavelengths in the LMC and SMC (Israel et al., 2010; Planck 
Collaboration et al., 2011a). Although, the COBE data of our 
Galaxy presented a submm excess (Reach et al., 1995), the inten- 
sity of this excess is much more pronounced in low-metallicity 
systems. Several explanations are in competition, for the origin 
of this excess: (/) very cold dust, in dense clumps, accounting 
for a large fraction of the dust budget of the galaxy (Galliano 
et al., 2003, 2005; Galametz et al., 2009; O'Halloran et al., 
2010); (ii) temperature dependent submm emissivity (Meny 
et al., 2007); rapidly spinning grains in addition to another 
component (Bot et al., 2010; Planck Collaboration et al., 201 la). 

The first Herschel observations of the LMC showed that the 
slope of the submm SED was flatter than in the Galaxy (Gordon 
et al., 2010). Meixner et al. (2010) showed that modefling this 
SED with standard Galactic grain properties required too much 
mass, and therefore concluded that it required modified grain 
optical properties. Roman-Duval et al. (2010) confirmed the 
Bernard et al. (2008) "far-IR excess" toward several molecular 
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clouds. Hony et al. (2010) demonstrated the complex structure of 
two massive starforming regions. Gordon et al. (2010) reported a 
SPIREsoo^m excess, which is likely the rise of the submm excess 
previously discussed. 

The unprecedented sensitivity and wavelength coverage of 
Herschel, at far-IR/submm wavelengths, allow us, for the first 
time, to study in detail processes that were previously glimpsed 
at. With a rigorous method, accounting for the different sources 
of error, it is now possible to unveil the systematic effects in- 
herent to SED modelling. The common assumptions, concern- 
ing the universality of dust properties, the accuracy of gas mass 
estimates, and the homogeneity of gas-to-dust mass ratios can 
be confronted by data. From a technical point of view, all the 
processes that have been previously described here define the 
required model parameters, as well as the unknowns when inter- 
preting the IR/submm emission of the LMC. 

For that purpose, the present paper is organized as follows. 
Sect. 2 presents the data set upon which we have based our 
analysis, and discusses the reference observational quantities we 
consider In Sect. 3, we present our SED model, and the way we 
rigorously propagate the various sources of observational errors 
through the entire fitting procedure, in order to provide reliable 
errors on the derived physical parameters. Sect. 4 attempts to rec- 
oncile different interpretations of the peculiar far-IR properties 
of the LMC: modified grain composition and/or undetected gas 
reservoir We end by a discussion on the origin of the SPIREsoo^,,,, 
excess. Sect. 5 synthesizes the paper and emphasizes the conse- 
quences of our findings. The various appendices give details on 
technical points that would otherwise alter the flow of the dis- 
cussion, if they were included in the main text. 

2. The Data Set 

The data set we are using are the science demonstration (SD) 
HerschellSPlKYi observations of the LMC (Meixner et al., 
2010), together with its Spitzer/IRAC and 5/7!teer/MIPS data 
(Meixner et al., 2006), and some ancillary data. These SD 
Herschel data cover only one fourth of the LMC. Although the 
complete PACS and SPIRE maps have now been obtained, we 
have performed our analysis on the sole SD strip, since we want 
to demonstrate general effects on the dust estimate, that do not 
require the totality of the LMC. 

2.1. Herschel Observations and Data Reduction 

We use the SPIRE maps presented by Meixner et al. (2010). We 
point the reader to this paper for the detailed description of the 
data reduction. The SD SPIRE data of the LMC cover a 2° x 8°, 
at 250, 350 and 500 yum. The extended source calibration was 
performed assuming SPIRE beam areas of 395, 740 and 1517"^, 
at 250, 350 and 500 fim, respectively. 

A background was subtracted, taking as a reference the two 
outer edges of the strip. Those edges are supposed to be out of 
the LMC. The same regions are considered for the background 
subtraction at other wavelengths, and for the gas maps. 

First, as discussed by Bernard et al. (2008), there are resid- 
ual foreground Galactic filamentary structures. We used the Hi 
map, whose velocity range corresponds to the Galaxy, in order 
to quantify the contribution of these fluctuations. Converting the 
Galactic H i column density into Galactic IR emission (using the 
Zubko et al., 2004, model), we find that the (non-subtracted) 
foreground accounts for ^ 15 % of the IR power of the strip. 
When this foreground is subtracted, the remaining fluctuations 



are on average ^ 1 % of the IR power Therefore, this contami- 
nation is smaller than our uncertainties. 

Second, we note that the larger scale Planck images (Planck 
Collaboration et al., 201 la) show that the edges of our maps are 
not completely beyond the LMC. The Planck images show that 
there is outer emission, which is colder than that of the rest of the 
LMC, in particular to the South of the strip. This oversubtraction 
may bias the dust temperatures, that we will derive in Sect. 4, 
toward hotter emission. However, as will be discussed in Sect. 4, 
our results are based on an excess of cold emission in the SPIRE 
bands. Therefore, this potential oversubtraction is conservative. 

We have compared the data of the SD paper, which were 
a single scan, with the more recent, accurately calibrated data, 
which also include a cross scan. The relative difference between 
the integrated SPIRE fluxes of the two sets is less than 2 %. 

2.2. Spitzer Data 

The Spitzer data include the four IRAC bands (3.6, 4.5, 5.8, 
8.0 yum) and the three MIPS bands (24, 70 and 160 //m). They 
have been presented by Meixner et al. (2006). 

2.3. Gas Tracers 

To compare the gas and dust contents in the LMC, we have com- 
pleted our data set with observations of the atomic and molecular 
gas phases. These maps were presented by Roman-Duval et al. 
(2010). 

We use the [Hi]2icm map observed by Kim et al. (2003). 
Their original beam size is 1.0'. The total atomic gas mass in the 
strip is M^^^ - 7.14x 10^ Mq. The noise at the original resolution 
iscrH,(1.0')^ 1.07Mopc-2. 

We use the '^CO(J=l— >0)2.6mm map observed by Fukui et al. 
(2008, with the NANTEN telescope). The beam size is 2.6'. 
We assume a constant CO-to-H2 conversion factor of Xco = 
7 X 10^° Hcm"^(Kkms"')"', based on the virial estimate of 
Fukui et al. (2008). It accounts for the variation of the Xco fac- 
tor with metallicity, due to less efficient shielding of CO by dust 
and to a lower intrinsic C and O abundance We emphasize that 
this value of Xco is larger than in the Milky Way. It accounts 
for discrepancies of the CO-to-H2 conversion factor in regions 
where CO is detected. However, it does not take into account 
potential regions where large envelopes of H2 could be present, 
but where the CO would be massively photodissociated, and 
therefore not detected in emission by ground based radio tele- 
scopes. With this conversion, the total molecular gas mass is 
M^as = 2.23 X 10^ Mq. The noise at the original resolution is 

6.05 Mq pc"^. The total uncertainty in Alf^^ is dominated by the 
uncertainty in Xco itself. We will discuss that point in Sect. 4.3. 

The gas masses above include the mass of Helium and heav- 
ier elements. The mean atomic weight used is: 

/^LMC = TT-*-^ = 1-34, (1) 

1 - Jq - ZlmC 

where the mass fractions of HeUum and heavy elements are 
Fo = 0.248 and Zlmc = 0.5 x Zq = 8.5 x 10"^ respectively 
(Grevesse & Sauval, 1998; Pagel, 2003). The fraction of molecu- 
lar gas in the strip is Mf^J(M^^\ + Mf^J ^ 24 %. It is higher than 
integrated over the entire LNIC (- 10 %; Bernard et al., 2008), 
since the strip includes a large number of molecular clouds. 



4 



Frederic Galliano et al.: Variations on the Dust Mass Estimate in the LMC 



Label 


Angular resolution 


Linear resolution (Zpix) 


Total 


Number of pixels 
Dust modelling Comparison to H i & H2 


Rl 


42" X 42" 


10 X 10 pc 


193 X 794 


89245 




R2 


56" X 56" 


14 X 14 pc 


144 X 535 


50367 




R3 


112" X 112" 


27 X 27 pc 


72 X 267 


12750 




R4 


3.7' X 3.7' 


54 X 54 pc 


36 X 133 


3278 


2442 


R5 


7.5' X 7.5' 


109 X 109 pc 


18 X 66 


852 


645 


R6 


14.9' X 14.9' 


217x217pc 


9x33 


225 


169 


R7 


29.9' X 29.9' 


434 X 434 pc 


4x 16 


64 


52 


R8 


r X 1° 


869 X 869 pc 


2x8 


16 


15 


R9 


2° X 2° 


1738 X 1738 pc 


1x4 


4 


4 


R18 


2.17° X 8.7° o 4.3° 


1890 X 7560 pc <^ 3780 pc 


1 X 1 


1 


1 



Table 1. Characteristics of the maps modelled in the present paper Each label (Rl, R2, etc.) corresponds to the full data set {Spitzer, 
Herschel and gas maps), but with different pixel sizes. The resolution RIO corresponds to the integrated strip; it is a rectangle. That 
is the reason why there are two dimensions. In the text, we will refer to the geometric mean of these two dimensions. The three 
columns on the right side of the table list the number of pixels. The first {''totar) is the total number of pixels in the map, at a given 
resolution. The second {''dust modelling") is the number of pixels used for the dust modelling. It is the total number of pixels in the 
maps minus the pixels which give bad fits at Rl. When building lower resolution maps, we exclude these bad pixels. Consequently, 
for a given waveband, the sum of the pixel fluxes is rigorously the same at all spatial resolutions. The third number {"comparison 
to Hi andH2") is the number of pixels used when comparing the dust and gas properties. Since the '^CO(J=l— >0)2.6inm map is not 
defined on the entire strip, it corresponds to a smaller number of pixels. 



2.4. Exploring the Effects of Spatial Resolution 

All our maps are regridded and reprojected on a common frame. 
They have been convolved with various kernels, in order to 
match the spatial resolution of MlPSieo^m (38"). The full pro- 
cess is described in Gordon et al. (2010). 

We aim at studying the systematic effects that would bias 
the dust mass estimates. The spatial resolution is one of these ef- 
fects. Indeed, SED models are highly non-linear, since the power 
emitted by a grain at thermal equilibrium with the radiation field 

(temperature Tdust) is proportional to T^^^ ^ ^du.st' where /3 - 2 
is the standard "emissivity index" (Eq. A.l). Therefore, the sum 
of the modelled dust masses of regions is likely to be differ- 
ent than the modelled dust mass of the sum of the emissions of 
these N regions. In order to study this effect, we will model sev- 
eral maps of the same region, but with different pixel sizes. 

Table 1 lists the different resolutions. The highest spatial res- 
olution (Rl, 42", 10 pc) is slightly larger than the largest beam 
size (MIPSi6{)/jm)- We then construct each set of maps {Spitzer, 
Herschel and gas) by summing the pixels in a 2 x 2 pixel win- 
dow. We repeat this process until we reach the size of the full 
integrated strip (RIO). The resolution of the combined gas maps 
is lower than the dust maps. We will therefore not study the spa- 
tial distribution of gas-to-dust mass ratios at resolutions higher 
than R4. From a technical point of view, the '^CO(J=l— >0)2.6mm 
map does not cover the entire strip. Consequently, we define a 
subsample of pixels where both the gas and dust maps are de- 
fined. 



3. A Phenomenological Dust SED Model 

3.1. Motivations 

We have developed a model aimed at accurately fitting the ob- 
served mid-IR to mm SEDs of various regions of the ISM of 
the LMC. At the spatial resolutions we consider here (> 10 pc), 
the SEDs will likely be the combination of several regions with 
different physical conditions - photodissociation regions (PDR), 
diffuse ISM, etc. In principle, we should perform a radiative 
transfer model, in order to determine the irradiation of each mass 
element of the ISM, then compute its spectrum, and transfer the 



IR radiation to the observer. However, we do not have the nec- 
essary information on the detailed matter and stellar 3D distri- 
butions, at these spatial scales, to constrain this type of model. 
Moreover, such an analysis is unnecessary in our case. Indeed, 
we are interested in the dust mass estimate. The mass is domi- 
nated by large grains, at thermal equilibrium with the radiation 
field. The spectrum of these grains depends only on the stellar 
power they absorb (or on their equilibrium temperature), and not 
on the details of the stellar spectrum and spatial distribution. 

This is not going to be true for grains which are out of 
thermal equilibrium with the radiation field (with typical radius 
fl < 10 run), especially PAHs. The spectrum of these grains de- 
pends on the hardness of the radiation field which determines 
the maximum temperature up to which the grain is fluctuating. 
The fact is that most of the emission of these grains arises at 
short wavelengths {A < 50 pm), and is not contaminating the 
far-IR-to-submm SED. 

Finally, the regions we are considering here are optically thin 
in the IR. Some compact sources show signs of absorption in 
the mid-IR (9.8 and 18 pm silicate features, 15.2 jum CO2 ice 
feature, etc.; Kemper et al., 2010; Hony et al., 2011), but the 
bulk of the grain emission, at longer wavelengths is unaffected. 

Considering the previously exposed arguments, we could 
derive reliable dust masses by simply fitting the observed 
SEDs with a combination of several modified black bodies. 
Nonetheless, we still choose to fit a combination of realistic dust 
models, even if the very small grain (VSG) and PAH spectra 
will not be perfectly accurate, due to the lack of constraints on 
the radiation field hardness. This approach is not providing sig- 
nificantly better mass estimates, but is providing more reliable 
estimates of the physical conditions of the hottest equilibrium 
grains. This is crucial for the interpretation, as will be demon- 
strated in Sects. 4.3 and 4.4. In particular, our approach allows 
us to avoid unphysical fits where a hot equilibrium component 
will be fit in place of the PAH emission. 

Our phenomenological dust SED model can be decomposed 
in two levels: 

1. the dust SED of a mass element of the ISM, which is con- 
trolled by the microscopic grain properties; 
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2. the synthesis of several mass elements, to account for the 
macroscopic variations of the illumination conditions. 

Sects. 3.2 and 3.3 details these two levels. This model was pre- 
viously used, in particular, by Galametz et al. (2009, 2010), 
O'Halloran et al. (2010), Cormier et al. (2010), Hony et al. 
(2010) and Meixner et al. (2010). 

3.2. Dust SED of a Mass Element of the ISM 

Let's consider the SED emitted by an element of mass of the 
ISM, where we can assume that the starlight intensity is uniform. 
For simplicity, we assume that the starlight intensity heating the 
grains has the spectral shape of the interstellar radiation field 
(ISRF) of the solar neighborhood (noted U^(A); Mathis et al., 
1983). We parametrize its integrated intensity by: 



U 



/jm 
Jo.0912 fi] 



2.2 X 10-5 Wm-2 



(2) 



The value U - I corresponds to the intensity of the solar neigh- 
borhood. In these conditions, large interstellar silicates have an 
equilibrium temperature of ^ 17.5 K. It is possible that the ISRF 
of the LMC differs from t/°, since it has younger stellar popu- 
lations. This spectral shape is also likely to vary spatially within 
the LMC, being harder in transparent regions, and redder in 
dense regions. However, as explained previously, this will im- 
pact only the PAH and VSG spectra, which are not determinant 
to estimate the total dust mass. 

In the same way, we are adopting the Galactic framework, by 
using the grain properties of Zubko et al. (2004). We choose the 
bare grain model, with solar abundance constraints (BARE-GR- 
S). The abundance and size distribution of each grain species 
was constrained by fitting the IR emission, UV-visible extinction 
and elemental depletions of the Galactic diffuse ISM. We have 
updated this model with the new Draine & Li (2007) PAH optical 
properties, that includes more accurate band profiles, based on 
Spitzer spectra. 

The optical properties and enthalpies considered here are 
summarized in Table 2. The grain cross-sections are computed 
using a Mie code, and following the method of Laor & Draine 
(1993, Sect. 2.2; App. A of the present paper). The temperature 
fluctuations are computed for each grain size of each compo- 
nent, and for each starlight intensity, using the transition matrix 
method (Guhathakurta & Draine, 1989). 

The specific monochromatic luminosity of a mass element 
of ISM exposed to the starlight intensity U is: 

/J"^'(f/, A) = /pAHC(f/> ^) + U.rhlT\U, A) + fsillfiU, A), (3) 

where /pah. /carb and fsn are the mass fractions of PAHs (charge 
fraction of 1/2), graphite and silicate (/pah + /carb + /sii = 1), 
and ll^'^iU,A), ll'''°{U,A), lf{U,A) are the corresponding size 
distribution integrated specific monochromatic luminosities. In 
this paper, we keep the mass fractions to the Galactic values 
(Zubko et al., 2004), except /pah that we vary to fit the observed 
IRACg^jm band. 

Throughout this paper, we will systematically compare the 
two following dust compositions, in order to study a possible 
evolution of composition between the Milky Way and the LMC. 

1. "77ie standard model" (hereafter labeled Std) is the origi- 
nal Zubko et al. (2004) Galactic grain composition, made 
of PAHs, graphite and silicates. The eff'ective submillimeter 



opacity of this dust corresponds to an emissivity index /? ^ 2 
(App. A). 

2. ''The AC modeV (hereafter labeled AC) is the ''stan- 
dard model", but replacing graphite by amorphous carbons 
(ACAR; Zubko et al., 1996). This substitution is arbitrary. 
The purpose of this model is to test realistic compositions 
having a higher effective submillimeter opacity {fi ^ 1.7; 
App. A), without violating the elemental abundances. 

We note that the heat capacities of amorphous carbons are un- 
known. We adopt those of graphite in replacement. It might be 
a crude approximation. However, this inconsistency will affect 
only the stochastically heated grains, which do not contribute 
significantly to the dust mass. 

The individual SEDs are shown in Fig. 1. The two models 
show similar features. 

1. The far-IR peak is dominated by grains in thermal equilib- 
rium with the radiation field. Their spectrum is roughly a 
modified black body. The peak wavelength of the emission 
shifts to shorter wavelengths when U rises, as the equilib- 
rium temperature rises. 

2. The mid-IR continuum, for U < lO'*, is dominated by small 
grains (radius a < 0.01 yum) and PAHs (prominent emission 
bands at 3.3, 6.2, 7.7, 8.6, 11.3 /^m), both out of equilib- 
rium with the radiation field. These grains are being heated 
by single photon events and the spectral shape of their emis- 
sion is independent of U. Their spectrum normalized to 
LiR (or U) is constant. The change of shape of the mid-IR 
spectrum with U is only due to the contribution of equi- 
librium grains at these wavelengths, when their temperature 
reaches Teq > 80 K (t/ > 10"*). In particular, the prominent 
9.7 jum silicate feature in emission dominates the mid-IR 
wavelengths, in this temperature regime. 

3.3. Synthetic Multi-Environment SED 

Each SED we model in this paper is likely to be the combination 
of the emission from regions with different physical conditions. 
To account for this diversity of conditions, we make the follow- 
ing assumptions. 

1. We assume that the dust properties are uniform within the 
modelled region: the size distribution, and mass fractions are 
constant. Only the starlight intensity varies. This is an ap- 
proximation. We will discuss in Sect. 4.3.6 potential local 
variations of the grain properties. 

2. The distribution of starlight intensities per unit dust mass, 
throughout the region, can be approximated by a power-law 
(Dale etal., 2001): 



dMn 



dU 



(4) 



This is an empirical prescription. Dale et al. (2001, Sect. 5.5) 
provide a physical justification of this formulation. However, 
its main advantage is that it allows for flexible parametrizing 
of the physical conditions. A more complex formulation is 
discussed in App. C.l. 



The total dust mass of each modelled region is therefore: 



-j 

Ju 



U^:„+AU 



dMd 



dU 



dU. 



(5) 
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Component Dielectric function Heat capacity Mass density 



Smoothed UV astronomical silicates 


Weingartner & Draine 


(2001) 


Draine & Anderson (1985) 


3.50 


X 


10^ 


kgm" 


-i 


Graphite 


Laor & Draine 


(1993) 


Dwek etal. (1997) 


2.26 


X 


103 


kgm" 


-3 


AGAR amorphous carbon 


Zubko et al. 


(1996) 


Dwek et al. (1997, graphite) 


1.85 


X 


103 


kgm" 


-3 


Neutral PAHs 


Draine & Li 


(2007) 


Dwek etal. (1997) 


2.24 


X 


103 


kgm" 


-3 


Ionized PAHs 


Draine & Li 


(2007) 


Dwek etal. (1997) 


2.24 


X 


103 


kgm" 


-3 



Table 2. Grain composition of our dust model. We give the references where the dielectric functions and heat capacities can be 
found, as well as the mass density of the material. 




Wavelength [/xm] Wavelength [/^m] 

Fig. 1. Uniformly illuminated dust SEDs, exposed to various radiation field intensities U. Each curve represents the sum of the 
emission by PAH*\ PAH^, carbon grains and silicates, exposed to U, in units of 2.2 x 10"^ Wm"^. The monochromatic luminosity 
is normalised by its integrated luminosity Lir. 



In addition, to subtract the stellar contribution from the mid- 
IR bands, we add a stellar continuum, parametrized by the stellar 
mass in the region M*: 



L*(A) ^M^x l*(A), 



(6) 



where I* (A) is the specific monochromatic luminosity of a 1 Gyr 
stellar population, synthesized with the model PEGASE (Fioc & 
Rocca-Volmerange, 1997). Since, this population is constrained 
mainly by the IRACs.e/jm and IRAC4.5JU,,, bands, the age of the 
populations do not have a big effect. On the other hand, this stel- 
lar mass is poorly determined and should not be trusted. The only 
purpose of this component is to give a better To accurately 
determine the mass of the stellar populations, we would need to 
take into account shorter wavelengths. This is not the purpose of 
this paper 

In summary, the total monochromatic luminosity of the 
model is: 



r mod 



(A) = 



l^^^^(U,A)x^^dU 



X It(A). 



(7) 



For each waveband /I,, where the observed monochromatic 
luminosity is L°''^(/i,), we compute the synthetic photometry, 
Uy°'^(Ai), by convolving the model with the instrumental spec- 
tral response, using the appropriate conventions provided by the 



user's manuals of each instrument. We minimize the x^, using 
a Levenberg-Marquart algorithm (Press et al., 1992). The;^^^ is 
weighted as follows: 



X 



{Lf\Ai) - L^°^{Adf 



(8) 



where AL^^(/l,) and AL^fiAi) are respectively the RMS and 
calibration errors of the waveband centered at wavelength Aj (see 
Sect. 3.4). We do not use the SPIREsoo^,,, flux as a constraint, be- 
cause of its excess relative to the model, as discussed by Gordon 
et al. (2010). Instead, we will study the behaviour of this excess 
in Sect. 4.4, in order to attempt to decipher its origin. 

Fig. 2 demonstrates this model on the integrated strip SED 
(RIO). The two panels highlight the degeneracy between the 
starlight intensity distribution and the submillimeter opacities, 
by comparing the fit of the two compositions to the same obser- 
vations. Having a flatter submillimeter opacity gives a fit with 
less massive, hotter dust populations. 

Table 3 summarizes the parameters of the model. In particu- 
lar, the last column of Table 3 describes the behaviour of the free 
parameters. When an SED fit is performed, a slight variation of 
one parameter value, due to observational errors, may system- 
atically be compensated by the variation of another parameter. 
Although the free parameters are rigorously independent, this 
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1 10 100 1 10 100 

Wavelength [yttm] Wavelength [/um] 



Fig. 2. Decomposition of the integrated strip SED into the individual uniformly illuminated SEDs. The grey circles and error bars are 
the integrated observed fluxes of RIO (Table 1). The total model (black line; Eq. 7) is the sum of the independent stellar component 
(grey filled area), and of the integral of uniformly illuminated dust SEDs (in colors; Eq. 4). There is linear gradation in U between 
colors for the uniformly illuminated SEDs. The starlight intensity U is in units of 2.2 x 10"^ Wm"^. The sum of these components 
is the black line. The green dots are the synthetic photometry (i.e. the model integrated in each instrumental filter). The left panel 
shows the "standard model", while the right panel shows the "AC model". To quantify the quality of the fits, the reduced chi square 
is^gj^ - 1.07 for the standard model" and;^^^ - 0.89 for the ''AC model" . 



effect may induce a correlation between these parameters. The 
Monte-Carlo error analysis that will be discussed in Sect. 3.4.2 
is a good way to quantify these correlations. The most striking 
example, in our case, is the correlation between the three param- 
eters controlling the starlight intensity distribution (a, and 
At/; Eq. 4). Moreover, these parameters do not have a physical 
meaning. Rather than discussing their values, it is more conve- 
nient to consider the first two moments of the starlight intensity 
distribution: 



1 



U X — dU 



which develop into: 



(U-{U}fx 



dMn 



dU 



dU, 



(9) 
(10) 



( 1-a (f/min + AU) 



u. 



2-a 



2-«(f/™„ + At/)'-'^-t/^;^ 

AU 

ln(C/mm+Af/)-lnf/n,i„ 

ln(U^in + AU) -InU^in 



if a 7t 1 & a 5t 2 



if a = 1 



c/„L-(f/™„ + Af/)- 



ifa = 2, 



and: 



cr'(f/) 



r3-a 



f l-a(t/,„„ + Af/)-"-f/4;„ 



' min 



3 - « (t/mi„ + At/)'-" - t/'T 



1 (U^in + AUf-Ul^^ 

[ 21n(C/n,in + A[/)-ln£/„ 



ifff^ 1. 



(12) 



We define the infrared luminosity as the power emitted by the 
dust: 



f 

Jo 



l;!"^'(v) dv. 



(13) 



Finally, the "gas-to-dust mass ratio" is the ratio of the total mass 
of gas (H I and H2, with Helium and heavy elements) to the total 
dust mass, in the same region: 



^dust 



(14) 



For comparison, the Galactic value is (Zubko et al., 2004): 
GI. - 158. (15) 

It is important to note that, for the discussion in Sect. 4, this 
value is consistent with the dust properties we use. Moreover, 
this model is consistent with the elemental depletion constraints. 



(11) 3.4. Rigorous Error Propagation 

Since our dust model is highly non-linear, it is crucial to rigor- 
ously propagate the observational errors through the entire fit- 
ting procedure, taking into account the fact that some errors are 
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Parameter 


Units & normalization 


Range of values 


Definition 


Comments 


1. Free parameters 


Mdust 


[Mo] 




[0,cx,[ 


Dust mass (all components; Eq. 5) 


Linearly correlated with ( U) 


/pah 


[4.6 %] 




[0,1/0.046] 


PAH-to-dust mass ratio (relative to 


Degenerate with the arbitrarily fixed 










Galactic) 


PAH charge fraction 


a 


... 




[1,2.5] 


Index of the power-law distribution of 


Correlated with U„ii„ and A 6' 










starlight intensities (Eq. 4) 






[2.2 X 10- 


-5 Wm-2] 


[10-2,10'] 


Lower cut-off of the power-law distribu- 


Correlated with a and At/ 










tion of starlight intensities (Eq. 4) 




At/ 


[2.2 X 10" 


-5 Wm-2] 


[1, 10' - U^in] 


Range of starlight intensities (Eq. 4) 


Correlated with a and Umin 




[Mo] 




[0,00 [ 


Mass of the old stellar population 


Poorly constrained; this parameter is 












used only to subtract the stellar contin- 












uum to the mid-IR bands; this compo- 












nent has no relation to the starlight in- 












tensity 


2. Derived parameters 




[Le] 




[0,cx,[ 


Infrared luminosity (total dust power; 


Almost model independent 










Eq. 13) 




m 


[2.2 X 10- 


-5 Wm-2] 


[10-2, 10'] 


Mass averaged starlight intensity 












(Eq. 11) 




cr(U) 


[2.2 X 10" 


-5 Wm--] 


[1,10' -<(/)] 


Second moment of the starlight intensity 












distribution (Eq. 12) 




Gdust 






[0,cx,[ 


Gas-to-dust mass ratio (Eq. 14) 


The Galactic value is G^^^^ - 158 












(Zubko et al., 2004) 



Table 3. Parameters of the model. The first category lists the independent free parameters. The second category lists the quanti- 
ties derived from these parameters. The last column ("comments") describes the behaviour of the parameters when fitting an SED. 
In absolute, M^^^t ™d {U} are independent. However, when fitting an observed SED, a slight variation of (U) will have system- 
atic consequences on Mdust, to compensate the variation and minimize the x^- The same comment applies to the other dependent 
parameters. 



independent and others are correlated. In this way, we will be 
able to quote consistent errors on the parameters. We first need 
to identify the various sources of error 

3.4.1. Sources of Observational Error 

For each wavelength, at each spatial resolution, we measure the 
noise of the map by taking the standard deviation of the pixel 
values in what is considered to be the background, i.e. the upper 
and lower ends of the strip. These values are given in Table 4. 
These errors are independent from one wavelength to the other, 
and from one pixel to the other. A simple check of the distribu- 
tion of pixel values shows that the uncertainty is well described 
by a Gaussian. 

The calibration error is the error on the flux conversion fac- 
tor. This error is therefore coiTelated between each pixel. It can 
be synthesized as follows. 

IRAC: the Icr calibration uncertainty is cr;,^! (IRAC) = 2 % 
(Reach et al., 2005). The different wavelengths are corre- 
lated. 

]VIIPS24;jm: Engelbracht et al. (2007) quote a Icr calibration er- 
ror of cr(,ai(MIPS24;,m) = 4 %. It is independent of the other 
wavebands. 

MlPSvo^m^ Gordon et al. (2007) quote a Icr calibration error of 

o"cai(MIPS7()^,„) = 5 % for the coarse scale mapping used for 

the SAGE observations. 
MIPSif,o;j/,,: Stansberry et al. (2007) report a Icr calibration eiTor 

of crcai(MIPSi60/im) = 12 %. This eiTor is coiTelated with the 

IVIIPSvop,,, error. 

SPIRE: Although Swinyard et al. (2010) report a calibration er- 
ror of crcai(SPIRE) = 15 %, it is necessary to decompose 
this error into its components (SPIRE consortium, 2010), as 



the SPIRE fluxes are the most crucial constraints on the dust 
mass and emissivity. 

1. The 3cr error on the calibration model is 3 x 
cr™'^'='(SPIRE) = 5 %, each waveband being correlated. 

2. The noise in the calibration observations (Ceres) are: 

- (r™|-(SPIRE250^,„) = 7 %; 

- cr™;-(SPIRE350;„„) = 12 %; 

- cr™;-(SPIRE5oo;.,„) = 6 %. 
These errors are independent. 

3. The error on the beam area are cr^^™(SPIRE) = 
-2 %/+5 % and they are independent. 

The latter SPIRE calibration errors correspond to the initial cal- 
ibration performed during the SD phase. We kept this error, in 
order to remain consistent with our data set. In any case, this es- 
timate is conservative, the new calibration errors being smaller. 

3.4.2. Monte-Carlo Iterations 

We propagate the observational eiTors detailed in Sect. 3.4.1, by 
performing Monte-Carlo iterations of each fit. More precisely, 
for each observed SED, we perform a large number (A^mc - 300) 
of fits of the SED with additional random perturbations. Fig. 3 
demonstrates the fits of the perturbed SEDs of the 4 pixels of 
the R9 map. Each model coiTesponds to one particular set of 
random perturbations. These perturbations take into account the 
two main sources of eiTors, as follows. 

1. The pixel noise at each wavelength (Sect. 4) is assumed to 
be a normal random independent variable. The noise is inde- 
pendent from one pixel to the other. 

2. The calibration error is assumed to be a normal random vari- 
able, with standard deviation and correlation between wave- 
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Filter [10'^ M Jy sr''^ 





Rl 


R2 


R3 


R4 


R5 


R6 


R7 


R8 


R9 


RIO 




19.8 


14.9 


7.4 


3.7 


1.86 


0.92 


0.40 


0.206 


0.104 


0.201 


IRAC4.5^„, 


2.61 


1.96 


0.98 


0.49 


0.244 


0.121 


0.053 


0.0271 


0.0137 


0.142 


IRACs.g^,,, 


3.8 


2.82 


1.41 


0.70 


0.35 


0.174 


0.076 


0.039 


0.0197 


0.239 


IRAC8^,„ 


3.7 


2.76 


1.38 


0.69 


0.34 


0.170 


0.074 


0.038 


0.0192 


0.45 


MIPS24;.,„ 


5.0 


3.7 


1.86 


0.93 


0.47 


0.230 


0.100 


0.052 


0.0260 


1.31 


MIPS70^,„ 


57 


42 


21.2 


10.6 


5.3 


2.62 


1.14 


0.59 


0.296 


16.0 


MIPSi60;j«! 


207 


155 


77 


39 


19.4 


9.6 


4.2 


2.14 


1.08 


19.8 


SPIRE250;jm 


248 


186 


93 


47 


23.3 


11.5 


5.0 


2.58 


1.30 


13.8 


SPIRE350^m 


180 


135 


68 


34 


16.9 


8.3 


3.6 


1.87 


0.94 


6.3 


SPIRE500;.m 


66 


49 


24.7 


12.3 


6.2 


3.05 


1.33 


0.68 


0.34 


2.73 



Table 4. RMS values for each filter and each spatial resolution. These values come from the standard deviation of the upper and 
lower borders of the maps. This RMS value represents the 1 x cr of a normally distributed random variable. 



lengths as described in Sect. 3.4. 1 . The calibration error from 
one pixel to the other is correlated. 

From a technical point of view, we generate the complete set 
of independent random variables necessary for the calibration 
errors, and keep them for our entire analysis. Indeed, one of the 
advantages of the Monte-Carlo technique is that it allows us to 
account for complex correlations between variables. As will be 
demonstrated in Sect. 4.2, the error on the ratio of two quanti- 
ties depending on the calibration error is often lower than the 
errors on each individual quantity. It is due to the fact that the 
calibration error cancels when considering a relative quantity. It 
is even possible to take this effect into account when comparing 
the results of two different models, as will be shown in Sect. 4. 1. 

Fig. 4 shows the statistical distributions of the main model 
parameters, corresponding to the fits of Fig. 3. The parameter 
value of each Monte-Carlo iteration is the sum of the parameter 
value of each pixel. The first panel of this figure demonstrates, 
in particular, that the dust mass has a ^ 50 % uncertainty, even 
for very high signal-to-noise ratio SEDs. In addition it clearly 
shows the asymmetry of the distribution, which is a result of 
the non-linearity of the model. The upper end of the statistical 
distribution of dust masses is less constrained than the lower end, 
because of the non-linear dependence of the dust mass with the 
temperature. 

Even with a fast running model, computing A^mc - 300 fits 
for each pixel is CPU intensive, and not necessary. Instead, we 
use an approximation to reconstruct the probability distributions 
of each parameter. This method is detailed in App. B. It consists 
of interpolating the pre-computed errors of 30 classes of SEDs, 
parametrized by their specific power (LiR/Mdust) and the RMS 
level. Fig. 4 compares this reconstruction method to the rigorous 
Monte-Carlo results. It succeeds in reproducing the correct cen- 
tral value and errors. In particular, it reproduces accurately the 
skewness of the probability distribution. 



3.4.3. Error display 

Throughout this paper, each time a numerical quantity is re- 
ported as X ^ a + b, the two quantities a and b will refer to 
the mean and standard deviation, or in other words: 



6i(^)i Qii^), Qi(^)^ corresponding to values of the repartition 
function of 1/4, 1/2 and 3/4, respectively: 



(16) 



On the other hand, for numerical quantities having a strongly 
asymmetric distribution, we will quote the error as X ^ a;^^. 
Only in this case will the quantities refer to the three quartiles. 



X - QiiX) 



-[Q2{X)-Qi{X)\ 



(17) 



The latter error (Eq. 17) corresponds to a confidence level 
of 50 %, by definition. It is also interesting to consider a higher 
confidence interval. Most of our error estimates are based on 
A^MC = 300 Monte-Carlo iterations. It would therefore be mean- 
ingless to go down to less than 1 / V^mc - 6 % error toler- 
ance. We therefore choose to quote the 90 % confidence level, 
defined by the range of the parameter values between 0.05 and 
0.95 of the repartition function. From a technical point of view, 
with A^MC - 300, the limits of this interval are simply the 15''' 
and 286''' ordered Monte-Carlo parameter values. We will note 
this interval: 



^ - Knf,^sup]90 %■ 



(18) 



On figures, the 50 % error bars will be displayed with a solid 
line, and the 90 % interval with a dashed line. 

We note that taking the median (Eq. 17) gives a central value 
very close to the best fit, while taking the mean (Eq. 16) gives 
a central value systematically shifted from the best fit. Fig. 4 
demonstrates this effect: the best fit value of the dust mass is 
8.62 X 10^ Mq, very close to the median of the Monte-Carlo 
iterations (8.76 x 10'' Mq). On the contrary, the mean, 1.03 x 
10^ Mo, is significantly higher. This is due to the skewness of 
the probability distribution of the parameter values. 

4. The Dust Mass Estimate and its Uncertainties 

4.1. Systematic Discrepancies Between the Two Models 

We have applied the model presented in Sect. 3, with the two 
compositions ("standard" and "AC"), to the maps listed in 
Table 1. The dust mass spatial distribution obtained with the 
two models are very similar. Fig. 5 shows the distribution of the 
pixel-to-pixel ratio of the dust masses {"standard" over "AC" 
models), and starlight intensities, at spatial resolution R4. These 
distributions are very tight. It shows, that the dust masses ob- 
tained with the "AC model" are systematically lower by a fac- 
tor of 0.38+IJ JJ^ than dust masses obtained with the "standard 
model", while the starlight intensities are systematically higher 
by a factor of 2.1;^y 2. These parameters are tied together. Grains 
with the "AC model" absorb more light than grains of the "stan- 
dard model" (App. A). They therefore require less mass to ac- 
count for the observed IR luminosity. 
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10 100 
Wavelength [/xm] 



10 100 

Wavelength [/Lim] 



Fig. 3. Demonstration of the Monte-Carlo method on the 4 pixels of the R9 map. The grey ellipses represent the unperturbed 
observations and their uncertainties: the filter band widths on the x direction; the ±lcr error bars on the y direction (dark grey); and 
the +3cr eiTor on the y direction (concentric light grey ellipse). The color lines show the A^mc - 300 model fits to the perturbed 
fluxes. We used the ''standard inodeF' for the demonstration. Some points (like MlPSvo^,,,) appears shifted from the models because 
of the color coiTection. The two shortest wavelengths of the [1,4] pixel are poorly fitted since the lRAC3,(,f,„i is only an upper limit. 
This might be the result of oversubtraction of point sources in this low surface brightness region. However, this discrepancy affects 
only the level of the independent stellar contribution. It does not aff'ect the longer wavelength fit. 



Thus, the results derived with the two compositions give sim- 
ilar trends, but the absolute value of their parameters systemati- 
cally differ. 



4.2. Bias Originating in tine Lack of Spatiai Resoiution 

In order to quantify the effect of the non-linearity of our SED 
model, we compute the total dust mass for each map listed in 
Table 1, with our two grain compositions. Fig. 6 shows the re- 



sulting trends of dust mass with spatial resolution. These trends 
are shown for both models. The mass at each resolution is the 
sum of the masses of all defined pixels. The top panels show 
the trend of the absolute value of the dust mass. The trends 
look similar for both models. They appear systematically shifted 
as discussed in Sect. 4.1. The error bar on each value is large 
and covers roughly the range of the trend. However, the rela- 
tive variation of the dust mass (bottom panels of Fig. 6) has sig- 
nificantly smaller eiTor bars. These trends are obtained by nor- 
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Fig. 4. Distribution of the main parameters of the perturbed SED fits of the fourR9 pixels (Fig. 3 ). For each parameter, we show the 
number of Monte-Carlo fits to the perturbed SED per bin of parameter value. The value of each parameter, for a given iteration is the 
combination of the values of the parameter for each of the 4 pixels. We quote both the mean and the median. The blue distributions 
correspond to the rigorous Monte-Carlo statistics, while the red distributions are from our reconstruction method (App. B). 



malizing each one of the A^mc Monte-Carlo results of a given 
spatial resolution, by the corresponding mass for the integrated 
strip (RIO). In that way, most of the calibration error cancels. 
The only remaining source of uncertainty is the intercalibration 
error between the various instruments and the RMS noise. The 
dynamics of the trend is unchanged. 

The relative trends of Fig. 6 show that measuring the dust 
mass of an integrated galaxy can give significantly lower values 
than performing fits of its spatially resolved regions, providing 
that the spatial resolution is fine enough. For the LMC, the spa- 
tially resolved dust mass estimate is ^ 50 % higher than the in- 



tegrated SED fit. This variation is not due to the noise, since the 
error bars are much smaller than the spread of the trend. Thus, 
the true mass is probably closer to the high spatial resolution 
value (Rl) than to the integrated flux value (RIO). It appears that 
there is a spatial scale where the trend stabilizes. This transition 
is an optimal resolution for our model, as it probably provides a 
coiTect dust mass, and the noise is lower than at Rl . This optimal 
spatial scale coiTesponds to R3-R4 (^ 27 - 54 pc; Table 1). 

The origin of this trend likely lies in the morphology of the 
ISM. It could be the result of the dilution of cold massive re- 
gions in hotter regions. It is possible that there is a typical spatial 
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Fig. 5. Comparison between the parameters obtained with the two models. The left panel shows the statistical distribution of the 
pixel-to-pixel dust mass ratio between the "AC model" and the "standard model" at spatial resolution R4. The right panel shows a 
similar ratio for the mass averaged starlight intensity {U). 



scale below which most of the cold regions dominate the SED 
of the pixels where they lie. It may correspond to the typical 
scale of molecular complexes. With only a few far-IR/submm 
constraints, it is difficult to account for these regions when mod- 
elling the integrated SED. 



Standard model 




10"^ 0.01 0.1 1 

Dust surface density E^u'^t [^o/pc^] 



Fig. 7. Pixel distribution of the dust mass surface density for two 
spatial resolutions. We plot only the "standard model". The fig- 
ure is qualitatively similar for the "AC model". The distributions 
are normalized. It shows that, at lower spatial resolution, the very 
high surface densities are missed (red filled area), and there is an 
excess of intermediate surface densities (blue filled area). 



To demonstrate the origin of this trend. Fig. 7 compares 
the pixel distribution of the dust mass surface density, at two 
spatial resolutions. Taking into account the statistical fluctua- 
tions, due to the low pixel number of R7, the two distribu- 
tions (Rl and R7) agree at low dust mass surface densities 
(^dust ~ ^-^^ pc"^)- However, at high dust mass surface den- 
sities (S^ll^j > 0.2 Mq pc"^), there is a higher fraction of pixels 
for the high resolution map. This fraction is compensated by an 
excess of intermediate surface density pixels at low spatial reso- 
lution (0.06 Mq pc-2 < I?^^^ < 0.2 Mq pc-^). Thus, the origin of 
the dust mass underestimation is due to the inability of our model 
to probe dense regions, at low spatial resolutions. At these spatial 
resolutions, the mass and average temperature of cold regions are 
biased by the contribution from hot regions present in the same 
beam, as the latter are more emissive. On the other hand, at high 
spatial resolution, cold regions tend to be better separated from 
hot regions and can therefore be modelled more accurately. 



4.3. The Gas-to-Dust Mass Ratio Crisis: Several Competitive 
Scenarios 

4.3.1. Metal Abundance Constraints 

To be rigorous, we first need to estimate the uncertainty on the 
metallicity of the LMC. Pagel (2003) compiles the literature for 
numerous elemental abundances in different regions of the LMC 
and of the Galaxy. We estimate the error on each element of 
Table 1 of Pagel (2003), by taking the dispersion of the different 
measures, excluding cepheids which are very dispersed and not 
available for all elements. These uncertainties are summarized 
in Table 5. 

We assume that we can reliably derive the metallicity by 
scaling the mass of O, C and N. In particular, this is supported 
by the fact that O and C are the major dust constituents. The 
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Fig. 6. Trend of the total dust mass with spatial resolution. The left panels show the trends for the "standard model", while the right 
panels show the trends for the "AC modeF. For each panel, the x-axis is the spatial resolution of the maps used to derive the dust 
mass. Each point of the trends corresponds to one of the maps listed in Table 1. For the two top panels, the y-axis is the total dust 
mass. This mass is the sum of the dust mass of each pixel. For each resolution, there are as many SED fits as the number of pixels 
listed in Table 1. The two bottom panels show the relative dust mass variation. It is normalized by the integrated strip (RIO). In that 
way, the calibration errors cancel, and the trend has smaller error bars. The dashed etTor bars display the 90 % confidence interval. 










C 




N 


12 + log (X///) Galaxy 


8.6- 


i.l 


8.4-8.6 


7.7 


-7.9 


LMC 


8.4- 


B.5 


7.1 - 8.1 


6.9 


-7.7 



Table 5. Literature compilation for the elemental abundances. 
The range of values reflects the dispersion of the observations 
listed in Table 1 of Pagel (2003). 



metallicity of the LMC we adopt is thus: 

16(0) +i2(C\ +i4(N) 



= OATSi- (19) 



This estimate is more accurate than simply scaling the oxygen 
abundance, as it is usually done. 

Let's check the physical consistency of our dust masses. The 
Galactic gas-to-dust mass ratio is G^^^^ - 158 (Zubko et al., 
2004). We emphasize that this value is consistent with the el- 
emental depletion patterns. Assuming that the dust-to-gas mass 
ratio scales with metal abundance {i.e. the dust-to-metal mass ra- 
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Fig. 8. Consistency of the gas-to-dust mass ratios. The two 
trends show the gas-to-dust mass ratio as a function of the spatial 
resolution, for each model. We also display the 90 % confidence 
interval for the standard model" , with dashed lines. The yellow 
and purple error bars, with the star symbol, represent the un- 
certainities on Gjj™; and G^^^^^, respectively (Table 5; solid line: 
50%; dashed line: 90%). The purple, green and yellow solid 
lines show the central values of G^^^^, G^^^^ and G^™;, respec- 
tively. The yellow filled area represents the gas-to-dust mass ra- 
tio range below the central value of the limit ratio Gj}";. 



tio is constant), the expected gas-to-dust mass ratio for the LMC 
is: 



^exp ggust ^ 339+55 



(20) 



Assuming that the mass fraction of gaseous heavy elements in 
the Galaxy is Zg ^ 0.017 (Grevesse & Sauval, 1998), the solar 
metallicity dust-to-metal mass ratio is: 



1 



^ 0.37. 



(21) 



It is therefore difficult to understand how the gas-to-dust mass 
ratio in the LMC could be lower than: 



^dust 



X ^ 125 



+20 



(22) 



without requiring a larger amount of metals locked-up in grains 
than what is available in the ISM. These values are summarized 
in Table 6. 



Quantity 



Value & Uncertainties 



■Zlmc/Zq 

^exp 
dust 



47+" "* 

"■^'-0.07 

125:^g 



[0.31,0.63190% 
[250, 503190% 

-dust —-15 [92,186190% 

Table 6. Reference metallicity and gas-to-dust mass ratios. 
These values are made explicit in Eqs. (19), (20) and (22). The 
convention for error display is defined in Eqs. (17) and (18). 



4.3.2. Preliminary: Global Analysis 



Resolution 



Rl 
R2 
R3 



R4 
R5 
R6 
R7 
R8 
R9 
RIO 



Gas-to-Dust Mass Ratio Gdust 
"Standard modeP' "AC model" 



88!i 



91 



[38,159190% 
[39,159190% 
[38,157190% 



"235^ 
237™ 



94 
107!'^ 
112+^^ 
110+1 
115! 



245!^^ 
250:| 
261+''^ 

299+11 
299+59 

-19 

020+40 



[86,410190% 
[86,410190% 
[84,410190% 



[78,400190% 
[78,400190% 
[81,410190% 
[85,420190% 
[90,440190% 
[88,430190% 
[91,430190% 



[39,158190% 
[40,160190% 
[42,162190% 
[47,174190% 
[49,179190% 
[47,180190% 

[50,173190% 

Table 7. Total gas-to-dust mass ratio, as a function of the spa- 
tial resolution, for the two models. This gas-to-dust mass ratio 
is the ratio of the total gas mass over the total dust mass, and 
not the averaged Gdust over the map. That is the reason why the 
high spatial resolutions are defined. The errors account for the 
uncertainties on both the dust and the gas masses. We give also 
the 90 % confidence interval. The convention for eiTor display is 
defined in Eqs. (17) and (18). 



Table 7 shows the gas-to-dust mass ratio at each spatial 
resolution for the two models. These ratios are displayed in 
Fig. 8. The total gas-to-dust mass ratios given by the "standard 



.+ L7 
'-LO 



model" are too low by a factor of G^^j^j/G^JJ^j(Rl) - 3.i 
[2.1,9.7]9o %. They are even lower than the physical limit by a 
factor of G;}™,/GS;,^,(R1) - 1.4!°;^ - [0.8, 3.6]9„ %. We empha- 
size here that there is no spatial coiTelation between the fore- 
ground Galactic Hi column density (Sect. 2.1) and the gas- 
to-dust mass ratio deficit of the "standard model". Therefore, 
the residual cirrus emission is not responsible for this deficit. 
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Statistically, the "standard modeV violates the elemental abun- 
dances with a probability of 80 %, while the "AC modeF is con- 
sistent. Therefore, if we assume that our gas mass is correct, then 
we can conclude that the properties of the "standard modeF do 
not apply to the LMC. This is one scenario. 

However, there is a second scenario: the discrepant gas-to- 
dust mass ratio obtained with the "standard modeF could re- 
sult from the underestimate of the total gas mass. In particular, 
the mass of molecular gas could have been underestimated. It 
is known that in low-metallicity environments, the H2 gas is 
not properly traced by '^CO(J=l— »0)2.6mm- In these environ- 
ments, the CO cores are thought to be much smaller relative 
to their H2 envelope (H2 being more efficiently self-shielded). 
This scenario is supported by the exceptionally high observed 
[Cii]i58;jm/'^CO(J=l— »0)2.6mm luminosity ratio in dwarf galax- 
ies (e.g. [Cn]i58;,m/''CO(J=1^0)2.6mm - 20000 in the LMC, 
compared to [C n]i58^ni/^^CO(J=l— >0)2.6mm - 4000 in normal 
metallicity galaxies; Poglitsch et al., 1995; Israel, 1997; Madden 
et al., 1997; Madden, 2000). In principle, the underestimation of 
the molecular gas mass, using the CO line, could be a factor of 
^ 10 - 100, in these environments (Madden et al., 2011). 

Let's assume that our dust properties are coiTect, and that 
^dutt ^^'^y '^^^ having underestimated the molecular 
gas mass. Noting Ch. the correction factor accounting for the 
hypothetical molecular gas not traced by CO and for the fact 
that our Xqo conversion factor (Sect. 2.3) might be wrong, the 
total gas-to-dust mass ratio would have to be: 

^Ch, - 10.1!^:^ = [4.0,30.6190%. 

In other words, to explain the discrepant gas-to-dust mass ratio 
of the "standard modeF, we would have to conclude that the 
molecular gas mass would have been globally underestimated 
by at least one order of magnitude. 

These two alternative scenarios are degenerate, when consid- 
ering only global values. We therefore need to take into account 
the redundancy provided by the spatial distribution of the gas- 
to-dust mass ratio, in order to sort these scenarios out. 

4.3.3. Spatial Distribution of tine Gas-to-Dust Mass Ratio 

Fig. 9 shows the map of IR luminosity for the "standard modeF. 
Since this quantity is the integration of the interpolated observed 
SED, it depends very little on the model; Lir distribution of the 
"standard modeF looks almost identical to the "AC modeF. That 
is the reason why we displayed it only for the "standard model". 
The upper panels of Fig. 10 show the spatial distribution of the 
gas-to-dust mass ratio with the two models. A few regions at 
the two ends of the strip exhibit noisy pixels. They correspond 
to very low surface densities. In general, there is no particular 
spatial correlation between the gas-to-dust mass ratio and the 
CO concentrations. The lower panels of Fig. 10 show the corre- 
sponding map of the mass averaged starlight intensity ([/), for 
both models. 

Fig. 1 1 shows the pixel-to-pixel distribution of Gdust for the 
highest spatial resolution where the gas maps are defined (R4). It 
appears that the distribution for the "standard model" is shifted 
to lower values, the pixels are systematically too low, com- 
pared to the expected value. In addition, we have built a his- 
togram of the pixels below an arbitrary column density (S^'^^^ < 
0.1 Mopc"^ and 1.^^^^ < 0.03 Mopc"^), in the two lower panels 



^iR [-^o ^] (standard model; R4) 2.75 




Fig. 9. Spatial distribution of the IR luminosity, for the "standard 
model", at resolution R4 (54 pc). The color image represents the 
IR luminosity map. The color scale is logarithmic. The white 
contours show the main CO concentrations from the Nanten map 
(Fukui et al., 2008). This contour level is chosen so that 90 % of 
the CO mass has a higher column density than this level. In other 
words, 90 % of the CO mass is in these concentrations. The map 
is almost rigorously identical with the "AC model". 



of Fig. 11. We have defined the surface densities at a given spa- 
tial resolution by dividing the mass in the pixel by the area of this 
pixel (Table 1): = ^^/'pix- lower panels of Fig. 11 

demonstrate that most of the pixels exhibiting high gas-to-dust 
mass ratios are located in regions of low dust column density, 
independently of the model used. 

Regarding this property, we need to study the variations of 
the gas-to-dust mass ratio as a function of surface density. 

4.3.4. Tine Correlation Between Gas and Dust Column 
Densities 

In this section, we analyze the correlation between gas and dust 
using a method similar to that of Planck Collaboration et al. 
(2011b). It is aimed at identifying the presence of "dark gas" 
(hereafter DG)', as a departure from the linear cotTelation be- 
tween gas and dust. We have defined bins of gas mass surface 
density (Ej^^ + S^Jj, so that the same number of pixels falls 
within each bin. We have then computed the average value and 
the scatter of the dust mass of the pixels within each bin. Fig. 12 
shows this binned trend on top of the pixel density plot. The 
binned dust mass surface density is fit with the following func- 

' In this paper, we choose the "dark gas" terminology, to remain 
more general. However, we acknowledge that CO-free molecular gas 
is not dark, since its dust grains are radiating and it likely emits copious 
amounts of [C ii] isg^m- 
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Cj„jt (standard model; R4) 3.454 Cj^jt (AC model; R4) 3.756 




Fig. 10. Spatial distribution of the main dust parameters, for the two models, at resolution R4 (54 pc). The color images of the 
two upper panels represent the gas-to-dust mass ratio map, for each model. The color images of the two lower panels represent 
the mass averaged starlight intensity map, for each model. The color scale is logarithmic. The white contours show the main CO 
concentrations from the Nanten map (Fukui et al., 2008). This contour level is chosen so that 90 % of the CO mass has a higher 
column density than this level. In other words, 90 % of the CO mass is in these concentrations. 

tion: for + C«'o2«i < and + C%I.^l, > 2^°, 

yHi I /^fit yH: _i_ yoff 

tit ^gas '-'CO^gas '^gas , , n , „ • 

- (/4j where the tree parameters are the following. 

^dust 
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Fig. 11. Pixel-to-pixel distribution of the gas-to-dust mass ratio for R4 (54 pc). The two histograms of the top panel are the pixel 
probability distribution of the gas-to-dust mass ratio for the two models. The spread of the distribution reflects the pixel to pixel 
spread fluctuations of the ratio. The actual error on the total Gdust is smaller than this spread (Table 7). For comparison, we show the 
values of G^^^j (Eq. 20), C^^^^ (Eq. 15), and Gjj™; (Eq. 22). In the two lower panels, each histogram is repeated and is decomposed 
according to its dust mass surface density. The hatched components correspond to lower surface densities (below the arbitrary limit 
indicated in each panel). The yellow and purple error bars, with the star symbol, represent the uncertainities on Gjj™; and G^^^j, 
respectively (Table 5; solid line: 50 %; dashed line: 90 %). 
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- GJj^^j is the fit gas-to-dust mass ratio. It assumes that the ac- 
tual gas-to-dust mass ratio is the same everywhere in the 
LMC. 

- C^Q is the fit correction factor of the assumed Xqo- Since 
we have assumed Xco - 7 x 10^" Hcm"^(Kkms"')"' 
(Sect. 2.3), then the value of the fit conversion factor is 
^co = C^o'^x lO^o Hcm-2(Kkms-i)-'. 

- accounts for a possible offset in the gas mass surface 
density compared to the dust mass surface density. 

- 1.^2 is the gas mass surface density above which the dark gas 
contributes. 

- is the gas mass surface density above which the molecu- 
lar gas is reliably traced by '^CO(J=l— >0)2.6mm, and the dark 
gas therefore does not contribute anymore. Above this value, 
the molecular phase is not "dark" anymore. 



Eq. (24) is not fit between 1.^^ and where the dark gas is 
assumed to contribute. The total dark gas mass is computed as 
the difference between the binned trend and the fit of Eq. (24). It 
is noted M°^. 

Fig. 12 shows the trend and the best fit for each model. 
Although the departure between the line fit and the binned trend, 
in the [Sg2^,Zg°] is visible, its deviation is lower than the typi- 
cal dispersion of the coiTelation. However, the parameter values 
and their uncertainties, given in Table 8, show this departure is 
significant. It appears that the gas-to-dust mass ratios derived 
from these simple fits are consistent with the pixel-to-pixel val- 
ues (Table 7). Therefore, this analysis tends to confirm that the 
''standard moder violates elemental abundances in the LMC, 
while the "AC model" is physically valid. The fact that the pa- 
rameter C[Iq is consistent with unity indicates that our adopted 
Xco conversion factor (Xqo - 1 x 10^" Hcm"^(Kkms ') '; 
Sect. 2.3) is probably not incorrect. Finally, the dark gas mass 
fraction, around 10 %, seems to be moderate. 

There are several hypotheses entering into Eq. (24). First, its 
interpretation implicitly relies on the translation between Sdust 
and Ay, the extinction magnitude in V band, because it de- 
fines the limits on the dark gas regime (2!?^*^ and Z^^; Eq. 24). 
However, at the spatial scales considered here (54 pc), molecu- 
lar clouds are not resolved, and several phases are mixed within 
each pixel. Therefore, the dust mass surface density is in princi- 
ple a biased estimator of the Ay. The value of Ay derived from 
the dust mass surface density, assuming a uniform dust distri- 
bution, is actually always going to be lower than the Ay of a 
molecular cloud that would lie in the pixel. Second, the approach 
of Eq. (24) is a perturbative method. It is correct only if the frac- 
tion of dark gas is small. In particular, this method would fail 
if dark gas was present at low Egas. Finally, it relies on the as- 
sumption of a uniform gas-to-dust mass ratio throughout the en- 
tire galaxy. In fact, our Gjust spatial distributions (Fig. 10) are 
not noise maps, they contain clear structures. Moreover, our er- 
ror analysis demonstrates that the amplitude of these structures 
is larger than the typical etTor bar on an individual pixel value. 
And the possible offset in background subtraction between the 
gas and dust maps (estimated in Table 8) is small enough to af- 
fect only the low surface brightness pixels. Consequently, this 
prompts us to further scrutinize the observed variations of Gdust 
as a function of the physical conditions. 




1 10 100 

Fig. 12. Correlation between the gas and dust mass column den- 
sities. Each panel corresponds to a model. The spatial resolution 
is R4 (54 pc). The colored area is the pixel density. The num- 
ber density of pixels for each [Xgas,Sdust] values is coded with 
the same color scale as in Fig. 13. The circles with error bars 
show the binned trend. The horizontal etTor bar displays the Egas 
bin width. The vertical error bar displays the dispersion in Edust 
within each Egas bin. The Egas bin sizes are chosen so that each 
bin contains the same number of pixels. The red line is the best 
fit of Eq. (24) to the binned trend. The vertical lines show the 
best fit values lP2 and E^S. The horizontal line shows the best 

fitvalueofofE°ft, = E°ff/GL. 



4.3.5. Variations of tine Gas-to-Dust Mass Ratio witli Physical 
Conditions 

In this section, we describe the observed trends of the gas-to-dust 
mass ratio as a function of several parameters. Fig. 13 shows the 
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Fig. 13. Pixel-to-pixel correlations between various tracers of the physical conditions and the gas-to-dust mass ratio. Results are 
shown for the two models, at spatial resolution R4. The color scale represents the density of pixels in various bins of the two 
parameters. The overplotted grey error bars are the trends binned over the x-axis parameter The bin size is defined so that the same 
number of pixels falls within each bin. That is the reason why the bins are very narrow in the center, and there are only a few points 
in the outer parts where the pixel density is low. The central value is the median gas-to-dust mass ratio, and the error bars account 
for the pixel-to-pixel scatter, which is similar to or larger than the typical intrinsic error bar of a single pixel. The yellow and purple 
dashed error bars, with the star symbol, represent the 90 % confidence uncertainities on Gjj™; and G^^^j, respectively (Table 5). This 
complex figure is explained in details in Sect. 4.3.5. 
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•2 [Mopc-2] 



[Mopc-2] 



CO 



"Standard model" 



131+3' 



1.57 



+0.45 
0.36 

6 9+' 

59+23 

1 42+0.41 



11.1 + 



-0.47 

;% 



"AC model" 



[55,2291,0% 


370 

-""-90 


[103,5901,0% 


[0.8,2.961,0% 


9 9T+0.34 
^•^-'-0.36 


[1.34,3.21,0% 


[3.5, 18.21,0% 


1 7 7+0.9 
i^.i j (J 


[9.6, 15.2],o% 


[53,981,0% 


54+3 


[53,981,0% 


[0.79,3.161,0% 


1 7Q+0.34 
^•^'-0.26 


[0.72,2.231,0% 


[0.45%, 24.5 %l,o% 


13.6+1^ % 


[3.8%, 25.9 %l,o% 



Table 8. Parameters of the fit ofLdust as a function of S^a, (Fig. J2; Eq. 24). The spatial resolution is R4 (54 pc). The uncertainties 
come from the Monte-Carlo analysis. The convention for error display is defined in Eqs. (17) and (18). 



pixel-to-pixel variations of the observed gas-to-dust mass ratio, 
as a function of several tracers of the physical conditions, for 
the two models. We note that the trends are similar with the two 
models, but the quantities proportional to the dust mass appear 
to be scaled down by a factor of ^ 0.38 with the "AC model". 
On the contrary, {U) appears to be scaled up by a factor of ^ 2.1 
(see also Fig. 5). 

The two left panels of Fig. 13 show the variation of Gdust 
as a function of the gas mass surface density Ec;a,s- There is a 
large scatter within each bin of Sgas, and the scatter is higher at 
low surface densities, since the signal-to-noise ratio is lower and 
offsets in the zero level of the gas and dust masses can impact 
the ratio. The gas-to-dust mass ratio at high gas mass column 
densities (Sgas > 10 Mopc"^) appears to deviate slighty from a 
constant value, as discussed in Sect. 4.3.4. However, this devi- 
ation is smaller than the dispersion. Overall, this trend does not 
show any general correlation between the gas-to-dust mass ratio 
and the gas mass surface density. The central values of Gdust are 
on average below the limit value with the "standard modeV\ but 
are found around the expected value with the "AC model". 

The two central panels of Fig. 13 show the variation of Gdust 
with the dust mass surface density Sdust- There is a clear anticor- 
relation between the two quantities. Although, the trend at low 
surface densities (Edust ^ 0.01 Mopc"^) is very scattered and 
sensitive to the offset 11°^^ (Table 8), the rest of the trend is clear. 
The fact that the gas-to-dust mass ratio coiTelates better with the 
dust mass column density than with the gas mass column density 
is the sign that the origin of the variations of Gdust is not directly 
linked with the observed gas content. Most of the trend is below 
the limit value, with the "standard model". On the contrary the 
values for the "AC model" are all realistic. 

Finally, the two right panels of Fig. 13 show the variation 
of Gdust with the mass averaged starlight intensity {U). There 
is a clear relation between the two quantities. Contrary to the 
central panels (Gdust vs. Sdust), these two quantities are indepen- 
dent. Moreover, (t/) is a specific quantity, while Edust is inte- 
grated over the line of sight. It appears that {U) is the parameter 
giving the best trend with Gdust, exhibiting larger dynamics and 
smaller scatter than the trends with Eg^s and Edust- We note that, 
with both models, the slope of the correlation (in bi-logarithmic 
representation) is not constant. Schematically, for (t/)^"* < 1 
(or {U)'^'^ < 2), Gdu.st {U). Then, for 1 < {Uf"^ < 4 (or 
2 < (U)^'^ < 8), the trend flattens. There is no significant vari- 
ation in this range. At larger values, Gdust (U) again, but the 
pixels are very dispersed in this range, and the statistics are lim- 
ited. 

In summary, the observed variations of Gdust do not appear to 
be significantly traced by the observed gas content. On the con- 
trary, these variations are well correlated with tracers of the dust 
physical conditions, in particular their irradiation conditions. It 
suggests that the origin of these variations is linked in some way 



with the radiation energy density. We are now going to explore 
the possible interpretations of these variations. 

4.3.6. Physical Interpretation: Dark Gas, Modified Submm 
Opacities or Enhanced Dust Condensation? 

To interpret the trends of Fig. 13, we need to list the possible 
physical processes responsible for variations of the observed 
gas-to-dust mass ratio, as well as the known biases affecting 
these trends. The known biases of our trends are the following. 

Bias 1. There is a potential offset between the gas and dust mass 
zero level, due to differences in background subtractions. 
This offset has an impact on low surface density pixels (or 
high ([/)). It is quantified in Table 8. 

Bias 2. There is a known bias of our SED fitting. At high av- 
erage starlight intensities ({U} > 10), our model tends to 
underestimate the dust mass. This effect is demonstrated in 
Fig. B.2 of App. B. An artificial correlation between {U} and 
Gdust is expected in this regime. 

The possible causes of variations of the observed gas-to-dust 
mass ratio are the following. 

Dark gas. As discussed in Sect. 4.3.2, low Gdust could be 
the sign of massive "dark gas" reservoir. As shown in 
Sect. 4.3.5, the fact that clear variations of Gdust are corre- 
lated with Edust, but not with Egas, suggests that the Hi and 
CO traced gas is blind to actual variations of Gdust- Assuming 
that our Xco is correct, the decrease of Gdust with Edust could 
be partly due to the presence of a dark gas component more 
massive than what was determined in Sect. 4.3.4. The latter 
could have been underestimated because it was lost in the 
scatter of the relation. 

Dust-to-metal ratio. Variations of Gdust could be due to local 
variations of the dust-to-metal mass ratio (Eq. 21). In prin- 
ciple, the gas-to-dust mass ratio could go down to Gjj™; 
(Eq. 22) and go up to infinity. Assuming an efficient mix- 
ing of freshly formed dust grains throughout the ISM, an in- 
crease of Gdust would be the sign of an enhanced destruction, 
mainly by SN blastwaves. This is probably the case at low 
surface density, although it is degenerate with biases 1 and 2. 
On the other hand, a decrease of Gdust would then be the sign 
of an increase of the condensation efficiency of metals onto 
dust. This effect is degenerate with the dark gas component. 
Assuming there is no dark gas, for this scenario to be real- 
istic, the metal condensation onto dust is expected to be en- 
hanced in regions where the gas density is higher. However, 
the left panels of Fig. 13 show that there is no clear coitc- 
lation between Gdust and the gas mass surface denstity. This 
effect should therefore be minor. Although the dust-to-metal 
mass ratio may vary in the LMC, our data indicate that this 
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potential variation is not at the origin of the trend, except 
maybe for low surface density. 
Emissivity. The variations of Gdust could be due to variations of 
the emissivity of the grains in different regions. This emissiv- 
ity variation would lead us to systematically misestimate the 
dust mass in regions where it would differ from the assumed 
opacity. These variations are poorly known and more diffi- 
cult to constrain. The most well-known of these processes, 
grain-grain coagulation, leads to an increase of the far-IR 
opacity (e.g. Stepnik et al., 2003). It happens preferentially 
in dense clouds. However, the left panels of Fig. 13 do not 
show a clear coiTelation between Gdust and Sga,. Therefore, 
this effect does not seem to be at the origin of the overall 
variations of Gdust- Another effect is the increase of the emis- 
sivity with grain temperature (e.g. Meny et al., 2007). This 
effect would lead us to overestimate the dust mass in regions 
where the dust is hot. It would therefore lead to an anticor- 
relation between the observed Gdust and {U}. Therefore, this 
effect can not be at the origin of our trend. Although this re- 
view is not complete, we will assume that local variations 
of the grain emissivity are not responsible for the observed 
variations of Gdust- 

Considering these previous effects, we now make several 
simple assumptions to interpret the trends of Gdust with {U} 
(right panels of Fig. 13). We assume that we can classify our 
pixels within the three following regimes, characterized by their 
average starlight intensity. 

A. We assume that there are embedded regions, characterized 
by low starlight intensities (1 < (t/)), where the observed 
gas-to-dust mass ratio is underestimated due to the presence 
of an undetected gas phase. The proportion of this over- 
looked gas rises, when (U) decreases. 

B. We assume that there are regions, characterized by moderate 
starlight intensities (1 < {U} < 10), coiTesponding to the Hi 
dominated regime, where there is no dark gas. The observed 
gas-to-dust mass ratio is assumed to be correct in these re- 
gions, and can be used as a reference. 

C. We assume that the gas-to-dust mass ratio may be higher in 
very diffuse regions, where shocks might have recently de- 
stroyed the grains. In these regions ({U) > 10), the biases 
of the modelling can also be significant, as demonstrated in 
App. B and Fig. B.2. An artificial correlation between {U} 
and Gdust is expected. 

Fig. 14 illustrates this scenario. The two panels show the 
binned trends of Fig. 13 with the starlight intensity, for the two 
models. The grey circles with error bars show the original trend, 
while the black circles with error bars show the trend corrected 
for the offset 1.°^^ derived in Table 8. We fit the latter. We decom- 
pose the trends into the three zones described above, by fitting 
the following function: 



log G 



dust 



log Gdust for <t^)d-k < (U) < {Uhis (25) 
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where Odark, Odiff, GJj'J^j, <t/>dai-k and <L')diff are free parameters, 
^dust "reference" gas-to-dust mass ratio, i.e. the gas-to- 

dust mass ratio in the diffuse ISM, and (f/)dark and {f/)diff are the 
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Fig. 14. Definition of the various zones for the two models, for 
R4 (54 pc). Each panel shows the binned trends of the right 
panels of Fig. 13 (grey circles with error bars). The fitted trend 
(black symbols with error bars) is the grey trend corrected for 
the offset 1.°^^ (Table 8). We define the three zones depending on 
their starlight intensities. Each zone corresponds to a different 
regime of the observed gas-to-dust mass ratio. The blue-green- 
red stripes are the envelopes of the fits of Eq. (25) to the trend. 
The yellow and purple dashed error bars, with the star symbol, 
represent the 90 % confidence uncertainities on Gjj"; and G^^^^, 
respectively (Table 5). This complex figure is explained in de- 
tails in Sect. 4.3.6. 



starlight intensities defining these zones (Fig. 14). Table 9 gives 
the main parameters of these fits for each model. 

This point of view leads to the same conclusion than in 
Sect. 4.3.4, namely that the "standard modeV is unphysical. 
Indeed, the derived reference gas-to-dust mass ratio, Gj' - 
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"Standard model" 


"AC model" 




Units 


^ref 
dust 


114+25 

1 QQ+O.Ol 
^•^'-0.21 
T, 4+0.6 


306+34 
3.6+"-5 
9 9+8:^ 


[2.2 X 10-5 
[2.2 X 10-5 


Wm-2] 
Wm-2] 



Table 9. Fitted parameters defining the three zones of Fig. 14. 
These values come from the best fit of the trends of Fig. 14 with 
Eq. (25), and the subsequent error propagation. 



4.3.8. Constraints on the Dark Gas Content 

Using the point of view of Fig. 14, we can put an upper limit 
on the dark gas content. Since we have shown in Sect. 4.3.6 that 
most of the variation of the gas-to-dust mass ratio could originate 
in an overlooked gas component, we can estimate the mass of 
this component, from this point of view. 



Mass fraction M^^, "'^ /(M^+ M^^,) 



114;^j4 (Table 9), in regions where dark gas is unlikely, is on 
average below the hard limit, Gjj™; - 117. On the contrary, the 
"AC model" gives a reference G^^^^^ - 306^3^ close to the ex- 
pected value, based on the metallicity of the LMC, G^^^^^ ^ 339. 
More importantly, this reference ratio does not violate the ele- 
mental abundances, with the "AC modeT. 

In summary, now that we have taken into account the various 
competing processes, we can safely conclude that the "standard 
moder, which works for the Milky Way, does not apply to the 
LMC. It means that the chemical composition of the grains in 
the LMC is systematically different than that of the Galaxy. The 
LMC grains have on average a larger submm opacity. Although 
the possible compositions are numerous, this work gives a plau- 
sible one, respecting the elemental abundances, based on optical 
properties of grains that have been observed in the laboratory 
(ACAR amorphous carbons in lieu of graphite; Zubko et al., 
1996). 



4.3.7. Remarks on the Submillimiter Emissivity Index 

First, as detailed in App. A, our "AC model" has a lower far-IR 
emissivity index (/3 ^ 1.7) than the "standard model" (J3 = 2). In 
our case, /? does not vary much with wavelength in the Herschel 
range. However, in general, y6 is a function of the wavelength, 
and there are grain species showing strong variations in this 
range (e.g. Quinten et al., 2002). 

Second, the emissivity of our grain model is what we call the 
intrinsic emissivity. On the other hand, the isothermal fit of an 
SED with free emissivity index provides the effective emissiv- 
ity index (noted jSeif. by us). The two are similar (J3 - ySeff.) only 
if the isothermal approximation is valid. On the contrary, if the 
spatial resolution is such that there is likely a mix of tempera- 
tures within the beam, then the two values are going to differ 
The SED will be broadened by the temperature distribution and 
ySetf. will likely be lower than /3. In the LMC, we are mostly in the 
second case. Due to this degeneracy, the shapes of the SED do 
not provide direct constraints on /3. Other constraints have to be 
invoked. In our case, we used the constraints on the gas-to-dust 
mass ratio to discriminate between different /?. 

We note that Aguirre et al. (2003), using COBE observations 
of the LMC, reached the same qualitative conclusion as derived 
here. They concluded that the grain opacities were different in 
the LMC than in the Milky Way, although they adopted a differ- 
ent set of grain species. However, these authors based their anal- 
ysis only on isothermal fits, and did not use the gas-to-dust mass 
ratio as a constraint. Therefore their demonstration was incom- 
plete. Similarly, Gordon et al. (2010, using Herschel data) and 
Planck Collaboration et al. (2011a, using Planck data) showed 
that the effective p of the LMC was around 1.5. However, they 
did not constrain the intrinsic /3. 

Thus, in light of these precisions, we have shown, for the 
first time, that the intrinsic far-IRJsubmm emissivity index of the 
grains in the LMC is lower than in the Galaxy. 



R4 
R5 
R6 
R7 
R8 
R9 
RIO 



"Standard model" 
87+52 % 
87+51 % 
85+49 % 
76!^ % 
69+1 % 
71+1 % 
64+f° % 



"AC model" 
52+4^ % 
50+49 % 
47+^« % 



42!| 
38_2 
37+ ■ 



% 
% 
% 



4i;| % 



Table 10. Properties of the overlooked gas component. This is 
the mass fraction of the optimistic estimate of the dark gas mass, 
defined in Eq. (26). This quantity is given as a function of spatial 
resolution for both models. 



To provide an optimistic estimate, we define the mass of 
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mass M^jjl. + M^^^, in order to match the expected gas-to-dust 
mass ratio G^^^'^ in regions where the uncorrected Gdust is lower 
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~,exp. 



(26) 



if Gdust > Gjy^j. 

Technically, this dark gas includes the CO-free H2, as well as a 
fraction of the CO-associated H2 where the Xqo factor is higher 
than our adopted value. This definition implies that this dark 
component covers most of the strip with the "standard model", 
while it is going to cover only a fraction of it with the "AC 
model". Table 10 gives the mass fraction of this component, for 
both models, at each spatial resolution. The most reliable value 
is probably the value at R4 of the "AC model". 

Therefore, comparing the conservative estimate of Table 8 
and the optimistic estimate of Table 10, with their respective un- 
certainties, we can bracket the dark gas mass fraction = 

M^S/iMf^^ + M^i,) ^ 10 - 100% in the strip, with the "AC 
model". 

For comparison, Bernard et al. (2008) found that this fraction 
was 182 %. The reason for this discrepancy might be due to the 
fact that the region we have modelled here is relatively poor in 
dark gas, compared to the entire LMC. 



4.4. What the SPIREsoo^,,, Excess Is Not 

As mentioned in Sect. 3.3, we have not fitted the SPIREsoo/jm 
excess present in our data. We have excluded this waveband, in 
order to avoid being biased by this effect, since its origin is still 
unknown. This excess can, a priori, affect our previous result 
concerning the grain properties, for two different reasons. First, 
this excess could originate in colder dust that our model has 
not accounted for Our dust mass would then have been biased. 
Second, if the origin of this excess is an unkown mechanism, 
not related to cold dust, it could also affect shorter wavelength 
bands, and bias our dust mass. 
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Although SPIRE5oo//m was not used as a constraint, for our 
dust SED model, it is possible to study the behaviour of its ex- 
cess with the physical conditions. 

4.4.1 . Systematic Analysis of tine Excess 
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Fig. 15. Different SEDs spanning the SPIREsoo/im excess range. 
These SEDs correspond to four pixels (at resolution R4; 54 pc), 
arbitrarily chosen. For each pixel the observations (open circle 
and eiTor bar), the model (solid line), and the model integrated in 
the broadband filters (filled circle) are displayed. For each SED, 
the corresponding dust mass surface density and the value of rsoo 
are given. 

Let's define the absolute and relative excesses of the 
SPIREsoo^m band, respectively by: 



A500 = vL°°'(SPlREsm^„,) - vLr''(SPIRE5oo;„J 
Lf XSPIREsoo^,,,) - C°d(SPIRE50(,;„„) 



'"500 



L°''^(SPIRE 



500//m 



(27) 
(28) 



Fig. 15 shows four SEDs with different values of the excess. 

Table 1 1 shows the value of the excesses, as a function of 
the spatial resolution. In general, the excess tends to be slightly 
smaller with the "AC model", since it has a slightly flatter submm 
slope. The relative excess tends to rise slightly when the spatial 
resolution decreases, but this trend is not statistically significant. 
We emphasize that rsoo being a relative quantity, most of the 
calibration error cancels, and the errors quoted in Table 1 1 are 
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7.6+H 



[-7,181,0% 
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[-7,161,0% 
[-7,151,0% 
[-4,131,0% 
[-3,111,0% 
[-3,111,0% 
[-1.9,91,0% 

Table 11. Absolute and relative SPIREsoo^m excesses as a func- 
tion of spatial resolution. This is the cumulative excess. The ob- 
servations are identical at all resolutions. Only the model and the 
error bars change slightly from one resolution to the other. The 
convention for error display is defined in Eqs. (17) and (18). 
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uncorrected between different spatial resolution. Most of the er- 
ror in this excess comes from the noise in the diffuse emission. 
The most reliable estimates of the excess are therefore those of 



RIO (integrated strip): 4*° = Ifi.l!* ^ % and r'^^^ 
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This excess is only marginally detected in the global SED. Since, 
the excess is relatively similar for both models, and that we have 
previously shown that the "AC modeV is more realistic, we will 
discuss the excess only for the "AC model", in the rest of this 
section. 

Fig. 16 compares the spatial distribution of the absolute 
SPIREsoo^m excess (Eq. 27, in colors), with the relative excess 
(Eq. 28, with contours). We first note that the absolute excess 
is not homogeneous, and contains structures associated with the 
known features of the LMC, in particular N44, and the north 
superbubble. We also note that the absolute value of the excess 
is much larger than the typical noise (Table 4). We emphasize 
that the map has been background subtracted in such a way that 
the fluxes on the north and south of the strip are zero. These ar- 
guments show that this excess can not be attributed to emission 
external to the LMC, like CMB fluctuations or cirrus foreground. 
The distribution of the relative excess rsoo also shows structures, 
that appear to be roughly anticorrelated with the IR luminosity. 

Let's look in more details at the correlation of the excess with 
different quantities, in order to decipher its origin. Fig. 17 shows 
the relation of the absolute excess to the gas phases. It shows that 
the absolute excess is not coiTelated with the molecular gas, but 
is coiTelated with the H i column density. It therefore suggests 
that the SPIREsoo/im excess is associated with the atomic medium 
but not with dense phases. To go further, we need to look into 
the various interpretations that have been proposed, for submm 
excesses in galaxies. 
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Fig. 17. Pixel-to-pixel SPIREsooiim absolute excess as a function of the column density of the atomic, and detected molecular phase. 
The spatial resolution is R4, the dust model is "AC. The color scale is identical to Fig. 13. The eiTor bars coiTespond to the 
binned trends. Eh^ is the column density of molecular gas derived from the CO line observations, but does not account for the dark 
component detected in Sect. 4.3. The yellow dashed line shows the zero value (no excess). 
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Fig. 16. Spatial distributions of the SPIREsoOf^m excesses. The 
spatial resolution is R4 (54 pc). The color image is the absolute 
excess Asoo, and the contours are the relative excess rjoo (the 
contours are labeled in percent). 



4.4.2. Consistency Test of Different Interpretations 

A submm emission excess has been reported in different sys- 
tems. We warn the reader that this excess extends up to cm 
wavelengths. It is still uncertain if the same physical process is 
responsible for the entire wavelength range, or if it is the combi- 
nation of several phenomena. 

Reach et al. (1995) first reported a long wavelength excess 
in the COBE/FIRAS spectrum of the Galaxy. This excess could 
be fit with a very cold dust component (4-7 K). However, the au- 
thors rejected this solution, since it was located at high Galactic 
latitudes. This excess was fit by Li & Draine (2001) invoking 
a long wavelength enhancement of the opacity of their silicate 
grains (see also App. A). Galliano et al. (2003) discovered a 
submm excess in the SCUBA and MAMBO observations of the 
blue compact galaxy NGC 1569. Then several studies report- 
ing similar excesses in galaxies, mainly low-metallicity dwarf 
galaxies were published, e.g.: Dumke et al. (2004), Galliano 
et al. (2005), Bendo et al. (2006), Galametz et al. (2009, 2010), 
O'Halloran et al. (2010). The relative intensity of this excess ap- 
pears to decrease with metallicity (e.g. the sequence SMC, LMC, 
MW). Planck Collaboration et al. (2011a) presented the global 
excess at long wavelength. 

The main interpretations of this excess, found in the litera- 
ture, are the following. 

1. Galliano et al. (2003, 2005) discussed the consequences of a 
very cold dust component, showing that it was not inconsis- 
tent if the emissivity index of these grains was /3 - I. This 
scenario would imply a small number (at most a few hun- 
dreds) of very dense, parsec-size clumps, containing a large 
fraction (40 - 80 %) of the ISM mass. 

2. Meny et al. (2007) showed that the excess was successfully 
fit with a physical model of temperature dependent grain 
emissivity. Their model implies that the emissivity index de- 
creases, when the temperature of the grain increases. 
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3. Bot et al. (2010) have shown that the long wavelength spec- 
trum of the excess renders the sole very cold hypothesis im- 
possible. It has to be combined with another process to be 
realistic: e.g. very cold dust and spinning dust grains (Draine 
& Lazarian, 1998). Alternatively, Planck Collaboration et al. 
(20 11 a) proposed that the excess in the SMC could be a com- 
bination of the Meny et al. (2007) grain model and of spin- 
ning grains. However, Planck Collaboration et al. (2011a) 
attributed the excess in the LMC to CMB fluctuations. Bot 
et al. (2010) explained the excess in the LMC and SMC with 
a combination of 12 K dust and spinning grains (Ysard & 
Verstraete, 2010). 

Here, we discuss rsoo in light of these scenarios. 

Since the actual value of this excess is very sensitive to the 
noise, let's zoom into the brightest star forming region of our ob- 
servations: N44 and its surroundings. Fig. 18 shows the central 
distribution of the excess around N 44, compared to the starlight 
intensity and the dust mass column density. Fig. 19 shows the 
coiTesponding pixel-to-pixel variations. It also demonstrates the 
typical error bars on individual pixels, showing that the excess 
is significant when looking at spatially resolved regions. Both 
(U) and Sdust are correlated with rsoo. However, the correlation 
is better with Sdust- This is confirmed by the spatial distribution: 
although, on average the excess corresponds to high starlight in- 
tensity regions, there is one spot north of N44, with no excess 
and high starlight intensities. Therefore, our observations indi- 
cate that the excess is primarily associated with diffuse regions. 
As a consequence, this excess can not come from contamination 
by the '2CO(J=5^4)520;,M line. 

In regard to these trends, the very cold dust hypothesis is 
very unlikely. Indeed, we would expect very dense clumps to be 
associated with dense regions. On the contrary, the observations 
show that most of these clumps would be present in the diffuse 
ISM, and very few or none of them in dense regions. However, 
this is only a qualitative argument showing that very cold dust is 
unlikely. To rule out this component, we develop the following, 
more quantitative point of view. In order to reach temperatures 
below 10 K, very cold dust has to be efficiently shielded. And, 
the dust responsible for the shielding is going to absorb the ambi- 
ent radiation and reach warm temperatures, emitting in the wave- 
length range that has been modelled in the previous sections. We 
can therefore estimate the mass of shielding dust necessary to al- 
low the very cold dust to reach low temperatures. By comparing 
the order of magnitude of the minimum mass of shielding dust 
to the observed mass in each pixel, App. D shows that very cold 
dust can reasonably be ruled out. 

Incidentally, by showing: (/) that very cold dust is un- 
likely to contribute significantly to SPIREjoo^m and ( ii) that the 
SPIREsoo^m excess is found mainly in diffuse regions, we have 
demonstrated that this excess does not have any impact on our 
main discussion about dust mass (Sect. 4.3). 

5. Summary and Conclusion 

In this paper, we have presented the modelling of the spatially 
resolved SpitzerflRAC, Spitzer/MIPS and Herschel/SPIRE data 
of a strip covering one quarter of the LMC. The purposes of this 
work was to: ( i) systematically study all the effects leading to in- 
accuracies or biases affecting the dust mass estimate of a galaxy, 
in the Herschel era; ( //) explore the peculiar ISM properties of 
the LMC. Our main results are the following. 

1. We have presented an empricial model to fit IR/submm 
SEDs. This model adopts realistic grain properties, and ac- 



counts for a possible distribution of starlight intensity in 
the region studied. We have described in detail the general 
propagation of the observational errors (noise and calibra- 
tion) through the SED fitting process. We have shown that, 
even when the signal-to-noise ratio is high, the eiTors on the 
dust mass are important and strongly asymmetric (typically 
~-25 However, we have shown that relative quantities ( ra- 
tio of two parameters) can have small error bars (typically 
-^7*%^) due to the partial cancelation of the correlated cali- 
bration errors. 

2. By modelling the same set of maps, but with different pixel 
sizes, we have shown that the lack of spatial resolution 
can lead to a systematic underestimate of the dust mass by 
^ 50 %. This bias could be the result of the veiling of cold 
components by the emission from warmer regions. Although 
the amplitude of this bias is probably specific to the type of 
environment found in the LMC and to the wavelength cover- 
age of our data set, we believe that this trend with spatial res- 
olution is general. Modelling the integrated SED of a galaxy 
leads to underestimating the dust mass. 

3. We have performed our SED fitting with two sets of grain 
composition, in order to explore the sensitivity of the dust 
mass estimate to the submm emissivity of the grains. We 
have shown that both compositions give equally good fits, 
but that the "standard" grain composition (graphite & sili- 
cate), that works for the Milky Way, fails to give physically 
realistic results in the LMC. It violates the elemental abun- 
dances. We discuss the fact that this discrepancy could ei- 
ther be due to modified submillimeter grain opacities (our 
AC model: amorphous carbon & silicate), or to the presence 
of an unaccounted for gas reservoir. In particular, we have 
shown that there is a degeneracy between the intrinsic grain 
emissivity and their temperature distribution, and that con- 
siderations on the mass are a powerful way to remove some 
of this degeneracy. 

4. The detailed spatial analysis of the variations of the observed 
gas-to-dust mass ratio has proven the need for grain opacities 
different than those of the Galaxy. More precisely, our analy- 
sis has demonstrated that grains in the LMC have on average 
a larger far-IR/submm opacity. We propose a physical dust 
model that is consistent with these properties (emissivity in- 
dex /3 - 1.7, and A:ab.s(160 pm) -1.6 m^kg"'), although we 
insist that this particular composition is not a unique solu- 
tion. 

5. We have shown that the mass averaged starlight intensity is a 
better tracer of regions where the observed gas-to-dust mass 
ratio is depleted. Comparing different methods, we have con- 
strained the amount of "dark gas" (i.e. unaccounted for by 
[Hi]2icm> and '^CO(J=l— »0)2.6mm)- This reservoir accounts 
for ^10- 100 % of the gas mass of the strip. Simply correct- 
ing the standard CO-to-H2 conversion factor is not sufficient 
to account for this dark gas reservoir, since it is abundant in 
regions where no CO detection is reported. 

6. We have analyzed the excess emission in the SPIREsoo/jm 
waveband, which has been reported before in the LMC and 
other low-metallicity galaxies. On average the amplitude 
of the relative excess is ^ 15 %, and it can vary spatially 
between and ^ 40 %. By scrutinizing its behaviour, we 
show that this excess is mainly antic or related with the den- 
sity of the ISM. We show that this excess can not be due to 
CMB fluctuations, nor to massive amounts of very cold dust. 
However, the nature of this excess remains unknown. 
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Fig. 18. Zoom on N44 and its surroundings. The color image is at R4, ■with the "AC modeV\ for both panels. The contours are 
(f/)^'^(/?4) for the left panel, and Z^^j(7?4) for the right panel. The white central star shows the location of the IR peak emission. 
The distributions are very similar, with the two models. 



Considering the uncertainties in grain optical properties, our 
study demonstrates that, when modelling the IR/submm SED 
of a galaxy, the derived dust mass can be considered as a con- 
straint rather than a result. Indeed, we have shown that standard 
Galactic grain properties were leading to unphysical masses. We 
therefore have selected another model which was giving reason- 
able results. However, going from this conclusion to providing 
the reader with an actual dust mass or gas-to-dust mass ratio for 
the LMC would be a circular process. It is now clear that grain 
properties can vary significantly from a galaxy to another. The 
range of gas-to-dust mass ratios that results from using different 
grain species makes the estimate of the absolute dust mass irrel- 
evant, without independent constraints. Such independent con- 
straints could be the elemental depletions, or spectroscopic in- 
formation on the precise chemical composition of the grains, or 
simultaneous extinction measurements. 



Appendix A: The Grain Properties of Our Models 

A.1. The Submillimetre Opacities 

The top panel of Fig. A.l compares the absorption opacities of 
the grain mixtures of our two models. These opacities are the 
sum of PAHs, carbon grains and silicates. The effective submil- 
limeter emissivity index /3 is defined by the logarithmic index of 
the opacity: 
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It is j0 ^ 2 for the "standard modeF and fi - l.l for the "AC 
model", although yS for amorphous carbons alone is even lower. 
Table A.l gives the submilimeter properties of our two models 
approximated by Eq. (A.l). 

The bottom panel compares the infrared emission of the 
compositions, for the same starlight intensity. It demonstrates 
that the "AC model" has more emissivity, especially in the 
submm. It therefore allows us to fit the same observed SED with 
less mass, and slightly hotter grains. 

Fitting the COB^/FIRAS high latitude Galactic cirrus, Li 
& Draine (2001) had to modify the imaginary part of the di- 
electric function of their silicate grains, at wavelengths greater 



Table A.l. Submillimeter properties of our dust compositions. 
The two parameters /3 and /Cabs(160 //m) are the parameters to 
approximate our dust opacities (Fig. A.l) with Eq. (A.l). 



than 250 i-im. This modification would have a very limited ef- 
fect on our conclusion, since our constraints go up to 350 jum 
only, where this increase is very limited (it is ±12 % in the 
250 yum < A < 1000 jum range). This modification actually 
lowers the emissivity between 250 and 850 fim, and increases 
it, at A > 850 fj.m. For our dust mixture, the amplitude of this 
modification is even lower, since the contribution of silicates to 
the far-IR is lower than for the Li & Draine (2001) model. The 
difference between the "standard model" and the one using the 
modified silicate of Li & Draine (2001) is invisible on Fig. A.l. 
The physical origin of this excess may be similar to the one we 
find here. However, its amplitude is much larger in the LMC, and 
manifests at shorter wavelengths {A> 100 fim). 

A.2. Size Distribution Considerations 

It is important to note that the size distribution used here does 
not include grains larger than 0.35 //m, and the contribution 
to the emission of grains larger than a > 0.1 fim is negligi- 
ble. Large grains tend to have a flat UV-visible opacity, and 
therefore, to have a lower equilibrium temperature than smaller 
equilibrum grains, exposed to the same ISRF. Therefore, adding 
larger grains would increase the submillimetre emissivity of the 
model, and allow us to fit the SPIRE fluxes, without having to go 
to low starlight intensities. However, having large grains or hav- 
ing colder dust would give the same discrepancies in terms of 
gas-to-dust mass ratios. Moreover, a significant excess of large 
grains would flatten the UV rise of the extinction curve, contra- 
dicting the observations of several lines of sight within the LMC 
(Gordon et al., 2003). 
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Fig. 19. Pixel-to-pixel variations of rsoo around N44, for the "AC model". This relation includes only the pixels of the window 
centered around N44 and displayed in Fig. 18. The spatial resolution is R4. The two top panels show the pixel density number 
coiTelations (colors; scale indentical to Fig. 13). The error bars represent the stacked trends in different bins. The central positions 
and the error bars are the median of the pixel distribution. They account for the scattering of the relation, but not for the eiTor 
on individual pixels. On the contrary, the two bottom panels show the same trends, but the ellipses are the median eiTor bars on 
individual pixels (dark: 50 %; light: 90 %). The excess is significant when looking at resolved scales. The yellow dashed line shows 
the zero value (no excess). 



Another feature of our model is that we have been forced 
to lower the abundance of non-PAH small grains (a < 10 nm; 
both carbon and silicate grains) by a factor of 2. Without this 
modification, we were not able to get good fits of the MIPS24;/«i 
of the diff'use regions. This modification has a very minor ef- 
fect on the dust mass (< 10 %), and it is systematic, thus it has 
no impact on our conclusions. However, this is puzzling since 
the fit of the extinction curves of the LMC indicates a larger 
fraction of small grains (e.g. Weingartner & Draine, 2001). The 



fact is that the Zubko et al. (2004, BARE-GR-S) model has a 
higher small grain contribution in MIPS24;jm than the Draine & 
Li (2007). On the other hand, replacing graphite by amorphous 
carbon, as shown in this paper, allows us to get rid of the 30 jum 
graphite feature, and decrease the contribution of small grains in 
the MIPS24^m band, making the 24 fim fit better without having 
to alter the size distribution. This is another indirect consistency 
check of the conclusion of this paper. 
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Fig. A.l. Comparison of the grain properties of our two models. 
The top panels show the opacity of the two mixtures. The grey 
dash lines show long wavelength fits of these opacities, with 
empirical laws oc i/. The bottom panel compares the SEDs of 
the two grain mixtures illuminated by the ISRF of the diffuse 
Galactic ISM (t/ = 1). 



A.3. Overview of the Derived PAH Properties 

Although the scope of our paper was not to discuss the PAH 
properties, their abundance is a natural output of our model. In 
this section, we summarize these results. 

Table A. 2 shows the PAH mass fractions for the two mod- 
els. This parameter is relatively well constrained and does not 
vary significantly with spatial resolution. Indeed, it depends es- 
sentially on the IRACs^m-to-total-lR luminosity ratio. The mass 
fraction for the "AC modeT is systematically higher, since the 
bulk of the dust is more emissive. 
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Table A.2. PAH mass fraction as a function of spatial resolu- 
tion. The quantity /pah (Eq. 3; Table 3) is the PAH-to-total-dust 
mass ratio, divided by the Galactic value (4.6 %). In other words, 
/pah = 1 in the Galaxy. The convention for error display is de- 
fined inEqs. (17) and (18). 
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Fig. A.2. PAH mass fraction, /pah- The map is shown for the 
"AC inodeF, at spatial resolution Rl (54 pc). 

Fig. A. 2 shows the spatial distribution of the mass fraction 
of PAHs, /pah (Eq. 3). We confirm the Paradis et al. (2009) re- 
sults showing an excess of PAH emission toward the stellar bar. 
However, /pah is biased by the fact that, in absence of detailed 
mid-lR spectrum, we have arbitrarily fixed the charge fraction 
to 1/2. The charge fraction controls the emissivity of the C- 
C modes (Galliano et al., 2008b). Therefore, it is difficult to 
uniquely interpret this excess emissivity by a local increase of 
the PAH abundance. 

Appendix B: Details Concerning tlie Error Analysis 

B.I. Classes of Observed SEDs 

The error analysis presented in Sect. 3.4 would imply having to 
fit the 156 577 pixels (Table 1) 300 times with the two models. 
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leading to a total of ^ 10^ fits. However, all these SEDs have a lot 
of similarities, and it is not necessary to perform the Monte Carlo 
iterations on each pixels. Here we describe the approximation 
method we demonstrated in Fig. 4. 



B.3. Biases of the SED Fits 
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Fig.B.l. Classes of observed SEDs for the "standard model". 
The SEDs are normalised by their integrated IR luminosity, Lir. 
The classes are ordered according to their specific power !Pdust- 
The classes derived with the "AC inodeT are very similar. Only 
the specific power is systematically scaled up by a factor of ^ 

i/o.38!;;:j;4(Fig.5). 



We start from the actual fit of the unperturbed SED of each 
pixel of Table 1. We order these SEDs according to their spe- 
cific emitted power !Pdust = i'lR/A^dust, making 30 logarithmic 
bins. Fig. B.l shows the classes of observed SEDs, for the "stan- 
dard model". At each waveband Aq, the central value is the 
median of the normalized monochromatic power A^)L,\^|Ll^ of 
each pixel within the considered bin of specific power. The error 
shows the dispersion of the pixels within the bin. On first ap- 
proximation, the specific power is proportional to the mass av- 
eraged starlight intensity {U). That is the reason why the SEDs 
of Fig. B.l are nicely ordered according to the wavelength peak 
of the dust emission. Notice that the dispersion at most wave- 
bands is small compared to the variation spread by the different 
classes. The only wavelength range where it is not true is the 
near-IR (/I < 5 ^va), where the stellar emission dominates. 

B.2. Class Interpolation and Error Estimate 

We perform the Monte Carlo error analysis of Sect. 3.4.2 on each 
SED class, for 6 different noise levels, spanning the whole range 
of observed signal-to-noise ratios, S/N. We perform this study 
with the two models. Then, for each pixel of Table 1, the errors 
on the parameters are determined by logarithmically interpolat- 
ing the precomputed errors in !Pdust and S/N. 

The validity of this interpolation method is demonstrated in 
Fig. 4. We show that this method reproduces well the central 
value and error bars of each parameter, since it reproduces well 
the skewness of the probability distribution. 
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Fig. B.2. Bias in the dust mass estimate as a function of the 
starlight intensity, for both models. In each panel, we plot the ra- 
tio of the median of the dust masses measured by fitting the ran- 
domly perturbed SEDs of Fig. B.l (M™f to the actual dust 

mass of the unperturbed SED (A^dust"'^' ^^^^ quantity is plotted 
as a function of the actual average starlight intensity {U). The 
color code of the classes is identical to Fig. B.l. We clearly see 
that for high starlight intensities, the dust mass is systematically 
underestimated by a factor up to 2. 



The large database of Monte-Carlo fits of the SEDs of 
Fig. B.l can be used to study the systematic effects of our 
method. Fig. B.2 shows the bias on the dust mass, as a function 
of the mean equilibrium grain temperature. 
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The two panels of this figure correspond to the two models. 
The jc-axis is the ''actual" {U}, i.e. the value of (U) derived from 
the unperturbed SEDs of Fig. B.l. Each value of {U} defines one 
of the 30 classes, since (U) oz ^P^usf The y-axes show the ratio 
between: ^™easured yyjjjcjj is fjjg median of the mass estimates 

dust ' 

derived from the fitting of the 300 Monte-Carlo perturbed SEDs 
of each class; and which is the dust mass coiTesponding 

to the unperturbed SED. In other words, quan- 
tifies the deviation from the dust mass derived from an SED fit, 
to its actual value, as a function of the starlight intensity. 

Fig. B.2 shows that up to (t/)^''' < 6 and {U}^'^ < 9, the 
dust masses derived from SED fits are not significantly biased 
(i.e. the deviation is smaller than the error bars). However, for 
starlight intensities higher than these values, the SED fit tends to 
systematically underestimate the dust masses (i.e. overestimate 
the gas-to-dust mass ratio). The amplitude of this effect can go 
up to a factor of ^ 2, for high (U). 

This is a demonstration of the effect invoked to explain 
the low gas-to-dust mass ratio, at high starlight intensities, in 
Fig. 14. Fortunately, this effect concerns only a regime contain- 
ing a small fraction of the pixels. Moreover, these pixels are the 
less massive, since they correspond to hot and diffuse regions. 
Therefore, this bias does not have a significant impact on our 
global dust mass estimate. 



Appendix C: Relevance of Our Starlight Intensity 
Distribution 
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Fig. C.l. Test of the stability of the starlight intensity distribution 
of Eq. (C.l). The parameter y is the mass fraction of the com- 
ponent in U^". The histogram is the distribution of the Monte- 
Carlo iterations of the fit of Eq. (C.l) to the integrated strip 
(RIO), with the grain composition of the ''standard model". The 
bimodality of the distribution is the sign of the instability of this 
formalism. 



C. 1 . The Unnecessariness of Adding a Diffuse Field 
Component 

We note that our approach of modelling IR SEDs with an em- 
pirical combination of starlight intensities, is common in the lit- 
erature. In particular, Draine et al. (2007) modelled the SEDs of 
the SINGS galaxies with an extra, uniformly illuminated, com- 
ponent: 



dU 



r^dust X 



(a -I) 



(£/n,i„ + Af/) 
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+ (\- r)Mdust X 6(U - t/mi„). 



where y is an extra parameter controlling the mass fraction of 
the U^" component. This formulation implicitly assumes that 
there are no massive quantities of dust colder than the diffuse 
ISM. One of the advantages of this extra component is to avoid 
the dust mass to diverge, in the absence of submm constraints. 
Indeed, as demonstrated by Galametz et al. (2011), modelling 
the SEDs of galaxies without submm data, using Eq. (4), leads 
to gross errors on the dust mass. However, in our case, we have 
the valuable submm constraints provided by Herschel. The rel- 
evance of this extra component can be tested using the Monte- 
Carlo method, described in Sect. 3.4.2. 

First, we tested the stability of this extra component. Fig. C.l 
shows the distribution of the Monte-Carlo values of y, when 
fitting the integrated SED strip (RIO), with Eq. (C.l) and the 
"standard" grain composition. This distribution is drastically bi- 
modal. Most of the perturbed SEDs are fit either without the 
extra 6(U - f/min) components (y ^ 1) - which is equivalent to 
our model - or with a sole uniformly illuminated SED (j ^ 0). 
In other words, the value of the parameter y is totally uncertain. 



Second, we tested the statistical relevance of adding this ex- 
tra component. To quantify this aspect, we performed a F-test. 
We computed the statistic (Bevington & Robinson, 2003): 
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where x\ii is the (Eq. 8) of the fit of a given perturbed SED 
of the integrated strip with the "standard model", and xlxtia is 
the corresponding value for Eq. (C.l); n - 10 is the number of 
wavebands, and mextra = 7 is the number of free parameters for 
the model of Eq. (C.l) plus the stellar component of Eq. (7). 

(C.l) 

Fy can be seen as a measure of how much the additional term 



has improved the value of the reduced x 



Statistically, F^ 



0.052+Qg27 ^ 0-67 is lower than the value of the F-distribution, 
with n - mextra -1=2 degrees of freedom, with a probability 



of exceeding F^ of 50 % (or Pf (0.67, 1, « - m^ 



1) = 50' 



Bevington & Robinson, 2003). In other words, adding the extra 
component of Eq. (C.l) does not significantly improve the;^'^. 
Our starlight intensity distribution (Eq. 4) is therefore statisti- 
cally more significant. 

Finally, we tested the conservativeness of the dust mass es- 
timate with Eq. (C.l). We have performed the fits of exactly 
the same perturbed SEDs that were used in Fig. C.l, with our 
"standard model" (Eq. 4). The ratio of the mass of our model 
to the mass obtained with the extra component is ■^diist/''^dust^ ~ 
l.OOO^QQ^y In other words, the two starlight intensity distribu- 
tions of Eqs. (4) and (C. 1) give statistically identical dust masses, 
within a few percents. 

In summary, for the analysis performed in this paper, the 
starlight intensity distribution of Eq. (C.l) is not physically mo- 
tivated, neither statistically stable nor relevant, and gives iden- 
tical results to Eq. (4). Our model has an appropriate balance 
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between free parameters and observational constraints. It is flex- 
ible enough, but it does not lead to overinterpreting the data. 



C.2. Comparison of Our Model with tine Isotinermai 
Approximation 

As justified in Sect. 3, our model accounts for a distribution of 
equilibrium dust temperatures within each pixel (Eq. 4). This 
formalism allows us to fit the submm slope more accurately than 
with a single temperature. As a consequence, our model pre- 
dicts more mass than a single black body fit having the same 
dust opacity. Since the single black body approach is still widely 
used, even in the Herschel era, we have performed a systematic 
comparison in order to quantify the biases of such an approach. 

To perform our comparison, we have fitted the R4 (54 pc) 
map with a single modified black body having the same grain 
opacity as our "standard modeV: /<-abs(160 jum) = 1.4 m^kg"', 
and /3 - 2 (Fig. A.l). For each pixel, the mass and the tempera- 
ture are free to vary, but /3 is kept fixed. We constrain this model 
with the MIPS70;,,,,, MIPSi60;.m, SPIRE250;.m and SPIREgoo;.,,. 
fluxes, weighted the same way as our complete model (Eq. 8), 
except that the weight of the MlPSvo^m fux is divided by 100. 
The purpose of this trick is the following: (i) when the peak 
of the SED is well constrained by MIPSi60/im, SPIREaso/^m and 
SPIRE350;jm (low temperatures), the MIPSvo^,,, flux has almost 
no impact on the x^- (H) when the peak is not well constrained 
(high temperatures), the MIPS7()^„, becomes important and tem- 
perature divergence is avoided. The mass for each pixel is noted 
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Fig. C.2. Comparison between our "standard model" and a sin- 
gle black body fit. The spatial resolution is R4 (54 pc). The x 
axis is the mass starlight averaged intensity {U} of the "standard 
modeT'. The y axis is the ratio between the dust derived with a 
sing le black body fit (m^^*^"^*) and with the "standard moder 

(^dust)' ^'^^ '■^^ same pixel. The color code of the pixel density 
and the binning of the trend are similar to Fig. 13. 



Fig. C.2 shows the ratio between the dust masses obtained 
with this isothermal fit and with our "standard model", as a func- 
tion of the mass averaged starlight intensity {U). It shows that for 
{U) > 1, the mass ratio is roughly constant between 60 and 80 %, 
but for {U) < 1, the ratio drops down to less than 30%. This 
simulation demonstrates that an isothermal fit will induce a bias 
in dust masses with average temperature. The regions with high 
starlight intensities {{U} > 1) are generally diffuse regions and 
their SEDs are almost isothermal. This explains why the mass 
ratio does not vary much with starlight intensity in this range. 
The average value of the ratio is around 70 % in these regions, 
which is reasonable considering the cmdeness of the approxi- 
mation. On the contrary, regions with low starlight intensities 
{{U} < 1) are generally dense regions. Since the mean free path 
of photons is much shorter in these regions, there is a significant 
mix of cold and hot components within each pixel. The isother- 
mal approximation is not valid anymore in this range of ([/), and 
the mass obtained with a single black body fit is biased. Namely, 
it drops by a factor of ^ 2 compared to its diffuse ISM value, 
down to ^ 30 % of the value of the "standard model". We obtain 
a similar trend comparing the "AC model" with a single black 
body fit having /(•abs(160 fj.m) = 1 .6 m^ kg"' , and yS = 1 .7, except 
that the values of {U) are systematically shifted by a factor of 
^ 2 (as demonstrated in Fig. 5). 

Finally, we note that the increase of dust mass at low starlight 
intensities compared to the isothermal approximation is unlikely 
a bias induced by the submm excess extending down to the 
SPIRE35()^„, band. Indeed, this excess is negligible in regions 
with low {U} (Sect. 4.4). It is prominent only in regions with 
high {U), where the trend of Fig. C.2 is not significant. 

Appendix D: Inconsistency of tlie Very Cold Dust 
Hypothesis 

As discussed by Galliano et al. (2003, 2005) and Galametz et al. 
(2009, 2010), the SPIREsoo^,,, excess could be attributed to very 
cold dust (VCD; T^^^ < 10 K). Simply considering the shape 
of the SED, this explanation is not unlikely. Such a component 
would produce a change in the submm slope of the SED. Indeed, 
warm and cold dust are heated by the transmitted stellar light, in 
different environments, at different optical depths. Their temper- 
ature distribution is therefore continuous. On the other hand, in 
order to reach very cold temperatures (Tdust ^10 K), the dust 
has to be fully shielded from stellar radiation. In these condi- 
tions, the dominant heating sources are the collisions and the IR 
radiation (galaxy and CMB) which do not depend on the opti- 
cal depth. It could reach temperatures lower than Tdust - 5 K 
(Galliano et al., 2003) only with difficulty. Very cold dust there- 
fore corresponds to a change of regime in the dust heating. If 
it exists, its temperature should not significantly vary with the 
optical depth in the cloud. Consequently its emission would be 
seen as a roughly isothermal component at 5 K < Tdust ^ 10 K. 
With sufficient mass, it would produce an excess emission, with 
a change of slope, at wavelengths 300 jum < A < 1 mm. 

Until now, we did not have the spatial resolution to explore 
this hypothesis. In this section, we do an order of magnitude 
estimate to test the likeliness of the VCD model. What follows 
is an update of the discussion at the end of Galliano et al. (2003). 
This is a complement to Sect. 4.4.2. 

Let's first assume that the excess originates in equilibrium 
grains having the UV-to-IR cross-section of silicates: 

V band opacity: /(-abs(/lv) - 130 m^ kg"' m n 

SPIRE5,)o^,„ opacity: ^abs(SPIRE5oo;,m)- 0.23 m2 kg-', ^ ' 
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and a submillimeter opacity index yS 



iVCD 



1 , and an equilibrium 



temperature of T^^^ ^ 10 K. Those are the most optimistic val- 
ues. If we can invalidate them, then the VCD hypothesis will be 
unrealistic for lower temperatures, steeper submm opacities, and 
larger UV cross-sections. 

Although the relative excess is stronger in diffuse regions, 
we find an excess in almost every pixel of the strip (Fig. 16). 
Let's assume that VCD lies in the core of very dense spheri- 
cal clumps in the diffuse ISM. Then, the ISRF they have to be 
shielded from is the typical ISRF of the diffuse ISM. For the 
more realistic model ("AC), this is basically {U) > 3. The opti- 
cal depth, assuming a slab extinction, to shield this dust, and to 
allow it to reach T^^^ is A^'^^ ^ 2.9. Since these clumps have 
to be small, and that we see the excess everywhere, it means that 
each clump is unresolved, even at Rl. Therefore, the diameter of 
these clumps has to be lower than the pixel size of Rl, which is 
/pix(/?l)- lOpc. 

This maximum size translates into a minimum density of 
n™°, in order to reach the required optical depth of A^'^^ ^ 2.9: 



4 VCD 

1.086 



^2.5x10'* Hem" 



(D.2) 



mil being the mass of an H atom. Notice that, since we are per- 
forming only an order of magnitude estimate, we use the slab 
extinction for a sphere. 

With {U} - 3, the temperature of the shielding dust, around 
the core, is TjJjJ.j''' ^ 21 K. Thus, a typical excess of r^QQ ^ 
15 % corresponds to a mass ratio, between the VCD core and 
the shielding ISM dust around, of: 



dust 



'"500 



/^"(SPIREsoo;,,,,) 



shield 



1 - ^500 '^"(SPIREgoO;,™) 



0.5, 



(D.3) 



Iv being the specific monochromatic power of the components. 
Therefore, for this picture to be correct, there should be at least 
^du'st''' of shielding ISM dust surface density, such that: 

This value is 2 orders of magnitude higher than the typical dust 
mass surface density, and one order of magnitude higher than 
the highest values (Fig. 13). Moreover, this argument is conser- 
vative, since the highest excesses are found in lowest surface 
density regions (Sect. 4.4). 

In summary, we have shown that the SPIREsoo/jm excess 
could not be accounted for by very cold dust, since it would re- 
quire a minimum mass of shielding dust, within each pixel, too 
large compared to the observed surface density. It is equivalent 
to note that, the excess being present in most pixels of Rl, we 
would need at least lO*" clumps of cold dust. It is more efficient to 
hide very cold dust in a small number of clumps. For example, 
Galliano et al. (2003) estimated that this number could be less 
than a few hundreds in the Magellanic dwarf galaxy NGC 1569. 
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